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Abstract: We present the two-loop QCD helicity amplitudes for quark-gluon scattering, 
and for quark-antiquark annihilation into two gluons. These amplitudes are relevant for 
next-to-next-to-leading order corrections to (polarized) jet production at hadron colliders. 
We give the results in the 't Hooft-Veltman and four-dimensional helicity (FDH) variants 
of dimensional regularization. The transition rules for converting the amplitudes between 
the different variants are much more intricate than for the previously discussed case of 
gluon-gluon scattering. Summing our two-loop expressions over helicities and colors, and 
converting to conventional dimensional regularization, gives results in complete agreement 
with those of Anastasiou, Glover, Oleari and Tejeda-Yeomans. We describe the amplitudes 
for 2^2 scattering in pure = 1 supersymmetric Yang-Mills theory, obtained from 
the QCD amplitudes by modifying the color representation and multiplicities, and verify 
supersymmetry Ward identities in the FDH scheme. 
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1. Introduction 



Recent years have seen rapid progress in our ability to compute two-loop matrix elements 
with more than a sing le kinematic variable |, |, |, |, |, |, |, |, |l|. The progress 
has relied in part on new developments in loop integration ||l^, |l^, 14 , |l^, |l^, and 



in understanding the infrared divergences of the theory |19|. The new two-loop amplitudes 
will be essential for reducing theoretical uncertainties in a number of physical quantities. 
(For a recent summary describing the various expected improvements see e.g. ref. pC|].) 
In particular, in jet physics an important source of theoretical uncertainty is from missing 

, which have been 



higher order corrections to the widely used NLO calculations [21, 22, 



crucial for matching theory to experiment. For the case of hadron colliders, one of the es- 
sential ingredients to performing next-to-next-to-leading order (NNLO) calculations are the 
parton distribution functions. Recently, global fits to the data for the parton distribution 



functions have been performed within an approximate NNLO framework |25]. There 
are, however sizable uncertainties associated with the experimental input to the parton 



distribution functions |26]. Nevertheless, an exact NNLO computation of jet production 
rates would be very welcome. Besides reducing the scale uncertainties for jet rates, such a 
computation will allow a better understanding of energy flow within jets, as a jet may con- 
sist of up to three partons at this order. For very large momentum transfer, improvements 



can also be made by resumming threshold logarithms |27]. 

Besides the two loop matrix elements, an NNLO calculation of two jet production 
at hadron colliders requires the tree amplitudes for six external partons |2^, 29 1 and the 
one- loop amplitudes for five external partons []30| , |3ll| , which have been known for some 
time now. Anastasiou, Glover, Oleari, and Tejeda-Yeomans have provided the NNLO 
interferences of the two-loop amplitudes with the tree amplitudes, for all QCD four-parton 
processes, summed over all external helicities and colors [Q, |5|, |6|. The helicity amplitudes 
for gg gg have also been computed, using the spinor helicity formalism [32|. In this 
paper, we present the qq gg and qg — > qg two-loop helicity amplitudes. The four-quark 
helicity amplitudes will be presented elsewhere. 

A useful property of helicity amplitudes is that they expose the full dependence on 
color and spin. Many formal properties of scattering amplitudes are simpler in a helic- 
ity basis and/or after color decomposition. Such properties include supersymmetry Ward 
identities |^, |4|, collinear limits ^ 36 1, and high-energy behavior |37, 38 1. The 
full color dependence is also useful for understanding the structure of the infrared diver- 
gences 111 



39|. 



This additional formal information is not necessary for the main phenomenological ap- 
plication, NNLO jet production in collisions of unpolarized hadrons. On the other hand, the 
helicity amplitudes can have phenomenological applications for jet production in collisions 
of polarized protons, as are being carried out at the relativistic heavy ion collider (RHIC) at 
Brookhaven. This may, for example, help to determine the poorly-known polarized gluon 
distribution in the proton |Eo|. Theoretical predictions of the relevant observables require 
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scattering amplitudes for polarized partons. Currently, predictions are available through 
NLO the helicity amplitudes presented here are a prerequisite for improving the pre- 
dictions to NNLO accuracy. Our results also serve as a check of the results of ref. ||5| for 
unpolarized qq gg and qg — > qg scattering. 

While preparing our results for publication, we became aware of a similar computation 



being completed simultaneously |^. Ref. |Q computes the one- and two-loop helicity 
amplitudes for qq — > gg and qg qg scattering, in the 't Hooft-Veltman (HV) scheme |45| 
(see below). A slightly different method from ours is used to extract the helicity ampli- 
tudes. We have compared our results for the two-loop amplitudes (and also the one-loop 
amplitudes), through the finite terms required for NNLO cross sections, and we are in 
complete agreement. 

In this paper we also describe the helicity amplitudes for gluino-gluon scattering in 
= 1 supersymmetric SU{N) gauge theory. Due to supersymmetry Ward identities p3| , 
these amplitudes are simply related to the = 1 gluon-gluon scattering amplitudes already 
presented in ref. Q. The supersymmetric amplitudes are a close cousin of QCD amplitudes, 
differing only in fermion multiplicities and non-abelian charge assignments: The A^ = 1 
supersymmetric amplitude are obtained from the QCD amplitudes by replacing the quarks 
with a single adjoint color representation fermion, effectively converting them to a gluino 
superpartner of the gluon. 

Because the scattering amplitudes possess both infrared and ultraviolet divergences, 
some care is required to ensure that the regularization procedure preserves supersymmetry. 
Several versions of dimensional regularization are commonly used for loop calculations in 
QCD. The widely used conventional dimensional regularization (CDR) scheme [43| breaks 
supersymmetry — it alters the balance between bosonic and fermionic states at order e, 
where e = (4 — Z))/2 and D is the number of dimensions. The CDR scheme is traditionally 
employed in calculations of amplitude interferences, such as in refs. |^, ^, |^, ^, In 
the helicity approach, the two commonly used schemes are the 't Hooft-Veltman (HV) 
scheme and the four-dimensional helicity (FDH) scheme p6| , |3^ . These schemes differ 
in the number of polarization states for unobserved gluons. The 't Hooft-Veltman (HV) 
scheme [45| contains 2 — 2e virtual gluon states (as does the CDR scheme), whereas the four- 
dimensional helicity (FDH) scheme [46, 34| assigns 2 states. The FDH scheme is related to 
dimensional reduction (DR) but is more compatible with the helicity method, because 
it allows two transverse dimensions in which to define helicity. It is also possible to define a 
scheme, labeled by a parameter 5ji, which interpolates between the HV {5ji = 1) and FDH 
{5f( = 0) schemes. We shall present the qq — > gg and qg qg amplitudes in this general 
6r scheme. A more detailed description of the differences between schemes, as well as a 
definition of the FDH scheme beyond one loop, has been given recently |34]. 

The supersymmetry preserving properties of the FDH scheme have been verified explic- 
itly at two loops, for particular helicity configurations of four-gluon amplitudes that vanish 
at tree level [||, §. In this paper, we explicitly verify that the FDH scheme preserves 
supersymmetry for the case of gluino-gluon scattering. Part of this check involves relat- 
ing gluino-gluon amplitudes to gluon-gluon amplitudes, for helicity configurations which 
are non-vanishing at tree-level. Such a test is somewhat more stringent than previous 
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ones P, 48, 34 1 because of the more intricate infrared divergences. The test additionally 
provides a nontrivial check on the calculation of the quark-gluon scattering amplitudes, as 
well as on the consistency of the FDH scheme. 

In general, scattering amplitudes in massless QCD possess strong infrared (soft and 
collinear) divergences. Using dimensional regularization, the amplitudes generically contain 
poles in e up to Catani has organized these divergences into a compact form predicting 
their structure [^]. We use Catani's formula and color space notation to organize the 
helicity amplitudes into singular terms (which do contain order terms in their series 
expansion in e), plus finite remainders. The precise form of the 1/e poles was not predicted 
a priori in ref. [^] for general processes at two loops. It is now clear, however, that these 
terms have a universal structure depending only on the external legs, based on explicit 
calculation ^, ^, ^, 0, ^, ^ |l^, |ll[ and on matching to resummed results [^] . 

For the quark-gluon scattering amplitudes discussed here, ref. previously computed 
the interference of the 1/e pole terms with the tree amplitude, summed over all colors 
and helicities. Here we extract the full color and helicity dependence of the 1/e pole 
terms. For the case of gg — > gg amplitudes, a term independent of color and helicity 
was found Q, plus a second "surprise" term with nontrivial color-dependence, which 
vanishes when the color-summed interference is performed. Here we confirm that a similar 
color-dependent term exists for the case of quark-gluon scattering amplitudes. We note 
that a term with similar color structure has been identified in contributions of one-loop 
factors for soft radiation to NNLO processes |49|. 

The conversion from one variant of dimensional regularization to another is well under- 
stood at one loop [46, 51]. At two loops, the ultraviolet shift in the coupling constant 
has been calculated for the commonly used variants of dimensional regularization p^ . 
However, the infrared aspects have not yet been fully understood for arbitrary processes. 
The scheme dependence of gg — > gg scattering amplitudes was studied in ref. There 
it was found that beginning at order 1/e^ the functions appearing in Catani's infrared 
decomposition are actually scheme-dependent. Here we present the scheme dependence 
(dependence on 6fi) for the qq — > gg amplitudes. The universal structure of the Catani 
formula for infrared divergences suggests that the conversion between schemes for any two- 
loop massless QCD amplitude would be controlled by the same set of functions that we 
uncover here and in ref. |p. However, the qq gg scheme-dependence we find here is 
significantly more intricate than that found for gg gg, and we have not yet identified 
the general pattern. 

The paper is organized as follows. In section ^ we review the infrared and color 
structure of one- and two-loop QCD amplitudes used to organize the amplitudes presented 
in this paper. In section ^ we describe the one-loop quark-gluon scattering amplitudes. 
Section 3.1 presents them in a form that is valid to all orders in e, in terms of integrals 
known through ©(e^). Their knowledge at this order is required for evaluating Catani's 
formula for the singular parts of the two-loop amplitudes through 0{e^). The finite one- 
loop remainder functions are given in section [3^ . The "square" of these functions, summed 
over colors and helicities, contributes to the NNLO cross section. In section |3^ , we describe 
how to carry out this sum, and compare the result to a similar sum performed in the CDR 
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scheme [|, ||. 

In section ^ we return to the two- loop ampHtudes. The finite two- loop remainder func- 
tions are presented in section 4.1 and appendix Some auxiliary functions for describing 
the shift in the finite remainder functions are given in appendix ^. The interference of 
the finite two-loop remainder functions with the tree amplitudes, summed over colors and 
helicities, also contributes to the NNLO cross section. In section 4.2 we describe the com- 
putation of this sum, and compare the results with those obtained in the CDR scheme in 
ref. §. 

In section ^ we discuss the = 1 super- Yang-Mills amplitudes obtained by modifying 
the QCD ones. The supersymmetry Ward identities are briefiy reviewed in section 5.1. In 



sections |5.3| and |5.4| we discuss the results for the one- and two-loop gluino-gluon scattering 
amplitudes for pure = 1 supersymmetric Yang-Mills theory, after first reviewing the 
infrared structure for the theory in section We verify that in the FDH scheme the 
two-loop amplitudes obey the expected supersymmetry Ward identities, which relate them 
to the gluon-gluon scattering amplitudes computed in ref. Q . In section ^ we present our 
conclusions. 



2. Review of infrared and color structure 

In this section we review the structure of the infrared singularities of dimensionally regu- 
larized one- and two- loop QCD amplitudes, using Catani's color space notation [|l^, as a 
prelude to presenting the finite remainders of the one- and two-loop qq gg and qg gq 
amplitudes. 

The two processes considered in this paper are 

Ai) + q{p2, A2) g{p3, A3) + 5((p4, A4) , (2.1) 
q{pi, Ai) + g{p2, A2) g{p3, A3) + q{p4„ A4) , (2.2) 

using a "standard" (not "all-outgoing") convention for the external momentum (pi) and 
helicity (Aj) labeling. The Mandelstam variables are s = {pi +^2)^1 t = {pi — ^4)^, and 
u = (pi -psf. 

We use dimensional regularization to handle both ultraviolet and infrared singularities. 
We consider a continuous set of schemes, labeled by a parameter characterizing the 
number of virtual gluon degrees of freedom circulating in loops. (Because we are computing 
helicity amplitudes, the number of external gluon states is fixed at two.) Specifically, when 
the trace of the Minkowski metric is encountered, we set 

ij% = Ds=i-2e6R, (2.3) 

corresponding to 2(1 — e6fi) gluon states in the loop. Setting 5r = 1 corresponds to the 
HV scheme, which is the most closely related to the CDR computation in ref. [||. Setting 
6r = corresponds to the FDH scheme, which has improved supersymmetry properties. 

The CDR and HV schemes imply the same coupling constant, the standard MS cou- 
pling, OsifJ-)- The coupling in a general 6r scheme is related to this coupling at NNLO 
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by 



6 ^ ^Mvr 



+ 



36 



12 



(1 - <5r) 



27r 



+ 



(2.4) 



Henceforth we will suppress the 6r index on ac,(/i). 

We work with ultraviolet renormalized amplitudes. The relation between the bare 
coupling a" and renormalized coupling as(/u), through two- loop order, is |18] 



2e 



asifJ') bo 
2tt e 



2lT 



(2.5) 



where /i is the renormalization scale, = exp[e(ln47r + V'(l))], and 7 = = 0.5772 . . . 

is Euler's constant. The first two coefficients appearing in the beta function for QCD, or 
more generally SU{N) gauge theory with flavors of massless fundamental representation 
quarks, are scheme-independent, 



lie A - 4TRNf 



l7Cl - {WCa + QCF)TRNf 



(2.6) 



6 ' 6 

where Ca = N,Cf = (iV^ - 1)/{2N), and Tr = 1/2. (Note that ref. [||] uses the notation 
/?o = 6o/(2vr),/3i = 6i/(2^)2.) 

The perturbative expansion of the qq gg amplitude is 



-"^-Mi!.U(/i;M) + 



(2.7) 



+ 



27r 



V 27r 



qq^gg ^ 
2 



-mSIpp(/";M) + o(«s(m)) 



where ■Mq^l^ggifJ-', {p}) is the L^^ loop contribution. The same type of expansion holds for 



the qg gq amplitude. Equation (|2.5D is equivalent to the following MS renormalization 
prescriptions at one and two loops, 



(1) = ^- 1 (Ij™ -^M 
-' ^qq-^gg ^qq~^gg ^ ■' ^ qq^gg ■ 



bo 



M 



(2) 

qq-*gg 



q-2 . .(2)unren 6o ^.-1 *.(l)unren ^1 A *.(0) 

-J^^qq^gg " -J^^qq^gg + I ^ ~ 26/ 



*gg • 



(2.8) 
(2.9) 



The infrared divergences of renormalized one- and two-loop n-point amplitudes are 
given by [|l|, 

|A^(i)(/x; M)),s = /(^He, M) \M'^\^^■. M))R.a + \M'^^^^{ii; M))aa , (2.10) 

|A^i')(/x;M))Rs = l^'\e,pi;{v}) |A^«(/.;M))aa (2.11) 
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where the "ket" notation \Aiij'\fi; {p}))n.s. indicates that the L-loop amphtude is treated 
as a vector in color space. The actual amplitude is extracted via 



A^„(ri,...,n"") = {ai,...,an\Mn{pi,---,Pn)) 



(2.12) 



where the aj are color indices. The subscript as. indicates that a quantity depends on the 

choice of regularization and renormalization scheme. The divergences of Adn^ are encoded 

in the color operator I^^^ , while those of A4n^ also involve the scheme-dependent operator 
r{2) 



In QCD, the operator I^^^ is given by 

1 „-ei>(l) " " 



2r(l -e) ^ 

^ ' i=i 



1 . JlI 

£2 ^ T? e 



2pi ■ pj 



(2.13) 



where Xij = +1 if i and j are both incoming or outgoing partons, and Xij = otherwise. 
The color charge Tj = {T^} is a vector with respect to the generator label a, and an 
SU{N) matrix with respect to the color indices of the outgoing parton i. For external 
gluons T". = so Tf = Ca = N, and 



75 



lie A - ^TnNf 
6 



For external fermions, the ratio 



7g_ ^ 3 

T? 2 



(2.14) 



(2.15) 



is independent of the representation. For quarks, Tf = Cp = {N^ — l)/(2iV); for gluinos 
Tf = Ca = N. 

The operator I^l is given by [|^] 

iS(£,m;{p}) = (^'"(t,f;{p)) + + 



+ 



r(i - e) 



+ i/[fl(e,/x;M), 
where the coefficient K^^ depends on and is given by [19, P 



(2.16) 



67_ 7r2 
18 6 



Ca - ^-^TnNf . 



(2.17) 



r(2) 



The function -Wrs. contains only single poles, and splits into two types of terms, 



4er(l -e) 



((2/f(2)+ 2/7(2)) 



1 + H 



(2) 



(2.18) 
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From the calculations performed here we find that the term proportional to the identity 

(2) (2) 

matrix in color space 1 contains the constants Hg and Hq , given by 

+ {-IcaCf + \cl + \cfTrN^ (1 - 5r), (2.19) 

+ ("i^l + l^A^R^f) (1 - ^«)- (2-20) 

This term survives the sum over colors, and the expressions for H^(f' and H^p in the 
HV scheme {5r = 1) agree, as expected, with previous color-summed results in the CDR 
scheme g |, |, |5|. The ijf ^ coefficient agrees for all 5r with that extracted from the 



99 ~^ 99 helicity amplitudes Q. 

The second term in ii"^^^ (e) has exactly the same type of nontrivial color and kinematic 
dependence found in the gg — > gg helicity amplitudes 1^], namely 



Ti • To , To • T, 



(2.21) 



with ln((— s)/(— t)) ^ Ins — ln(— t) — ivr in the s-channel, etc. (It might seem that the 
operator Ti • T2 is somewhat ambiguous for the qq — > gg process ( |2.lD , since the qq and gg 
pairs have different color quantum numbers. However, the difference between the qq and 
gg Ti • T2 operators is proportional to the identity — it is a difference of Casimirs — so 
the commutator is unambiguous.) 

In ref. it was observed for the gg — > gg amplitudes that the second term in eq. ( 2.1^ ) 



is independent of the helicity configuration, and is a nontrivial commutator matrix in color 
space. (The possibility of nontrivial color structure in H^'^\e) was pointed out in ref. [|ig|].) 

- (2) 

Now we can see that H is also independent of whether the external lines are quarks or 
gluons, which buttresses the suggestion P| that it is related to soft, not collinear, virtual 
contributions. A similar color structure emerges in a general analysis of the contributions 



of one-loop factors for soft radiation to NNLO processes [49|. Because of the commutator 



- (2) 

structure, H vanishes when sandwiched between tree amplitudes, after performing a sum 
over colors; hence it drops out of the color-summed interference of the two- loop amplitudes 
with the tree amplitudes [^. 

To proceed further, we wish to introduce an explicit color basis for the amplitudes, 
and also remove certain overall spinor product factors. To do the latter, we take the set of 
independent helicity configurations h to be 



h= 1 
h = 2 
h = 3 



q{Pi, +) + q{P2, -) ^ 5(P3, +) + 9{Pi, +) , (2.22) 
q{Pi. +) + (i{P2, -) ^ 5(P3, -) + 9{Pi, +) , (2.23) 
q{Pu +) + 9(P2, -) ^ 5(P3, +) + 9(P4, +) , (2.24) 
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Amplitude 


Value of h 


Permutation 


On — > Qo^ qT 


h = l 




On — > qT 


h = 2 




ot o<7 — > ot o7 


h = 2 


it <— > u\ 

L J 


9]^ C ^ 5^ 9a 


h = 1 


{t ^ u} 


(it at at it 


h = 5 




it a2 at it 


h = 3 




it at as it 


h = 3 




it a2 as it 


h = 4 





Table 1: Relations between a general helicity configuration containing q{pi, 
presented in the text, eqs. (2.22)-( 2.26| ). 



-) and the five 



h = 4: q{pi,+)+gip2,-) 
h = 5: q{pi,+) + g{p2,+) 



a{P3,-) + qiP4,+) , 
a{P3,+) + Q{p4,+)- 



(2.25) 
(2.26) 



Other configurations are simply related to these by symmetries. For example, the q{pi, —) 
amplitudes are obtained by parity (P); the qg gq amplitudes are related to qg gq by 
charge conjugation (C); q{pi,+) + q{p2,—) — > giP3,—)+giPi,—) is related to process ( 2.22| ) 
by CP; and q{pi, +) + g{p2, +) g{P3, —) + qip^^ +) is related to process ( 2.24 ) by time 
reversal (T). Maintaining helicity conservation on the quark line removes half of the con- 
figurations. Applying parity to the remaining configurations, in order to let the helicity 
of q{pi) be positive, there are four configurations each for qq —>■ gg and qg —>■ gq. Table || 
relates these eight configurations to the five represented by eqs. ( 2.22 )-( 2.26D . 

Helicity-dependent, phase-containing factors arise because we evaluate the amplitudes 



in the spinor helicity formalism |32]. For h = 2,4,5, we remove a factor Sh related to the 
tree amplitude. For h = 1,3, the tree amplitude vanishes, and we remove a factor related 
to the one-loop amplitude. We define 



Ss 



. (13) [34] 
^(2 3) (34) 
. (13) [3 2] 
'(2 3) (34) 



So 



(13)^ (2 3) 



(12) (2 3) (3 4) (41) 



5*4 



(13)'^ (4 3) 



(12) (2 3) (34) (41) 



{12f (4 2) 



(12) (2 3) (34) (4i; 



.(2.27) 



The spinor inner products |32, 29| are (ij) = and [ij] = (i+|j~), where \i ) are 

massless Weyl spinors of momentum ki, labeled with the sign of the helicity. They are 
anti-symmetric, with norm | (ij) | = | [ij] \ = 
prefactors enter polarized cross sections, 



Sij, where Sij = 2ki ■ kj. The squares of the 



1^1 1 



\Ss\' 



U IC |2 _ ^ 

U IQ |2 _ 

~ ) I '-'4 1 — ~^ 

S st^ 
where we have included a minus sign for all the qg 
cross sections. 



su 



(2.28) 



gq cases to ensure positive tree-level 
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The color decomposition of the amphtudes (or their finite parts, according to eqs. ( 2.1C| ) 
and ( |2.f 1| )) reads 



TW^^^ = Sh><Yl '^^''^ ^ ^'[^1 , /i = 1, 2, 3, 4, 5, (2.29) 



c=l 

where 



Trll] = (2^a3^a4)i2^ ^ rp^[2] ^ ^r^a4rj.a3y2^ ^ ^^.[3] ^ ^a3a4 ^2^ ^ /j = 1, 2, (2.30) 

and 

= (2-^32-^2^4^ ^ XrP] = (T«22-»3)i4^ ^ xr[3] = ^'^^^s ^ /i = 3, 4, 5.(2.31) 

Here are SU{N) generators in the fundamental representation, normalized according 
to the convention typically used in helicity amplitude calculations, tT:(T"'T^) = 6°'^. (The 
T" used in this color decomposition should not be confused with the appearing in 
whose representation depends on the external line; nor should they be confused with the 
generators for the quark representation, which have the more "standard" normalization, 
Tr = 1/2, as mentioned above.) The amplitude components depend on the Mandelstam 
variables, M^^^'''^ = Ml^^'^''\s,t,u). Often we will suppress the kinematic arguments, 
unless a permutation of them is involved. 

Note that we use a charge-conjugated fundamental index for quarks and anti-quarks 
in the initial state, i.e. ii for qi, and 12 for ^2- Also, amplitudes related to h = 1,2,3,4,5 
by charge conjugation (C) will naturally have the fundamental indices charge conjugated. 

In the basis ( |2.30| ) for qq — > gg, the matrix i"*-^-* is 

^^^(^) = -f(r3^^ (^-^^^ 
-2^c. + fe,)s + ive,sT Cg,(T-u) 

£95(3 -U) ^ggiS-T) im^g + NQS^ 

where 

-(^r -(^r -(^)*' 

e. = ^ + |, ^9 = ^ + ^' ^19 = 1(^1 + ^^ (2-34) 

V = N"^ -1. (2.35) 

The corresponding operator for qg gq in the basis ( p. 31 ) is obtained by exchanging S 
and T in eq. ( 2.32| ). 

The tree amplitudes in the color basis ( |2.3C1| ), ( |2.31 ) are given by M^^'^^\ where 
= '[^1 = 0, c = l,2,3. 
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»^(0),[3] 


= 0, 






n = 1,2, 6, 4, 5, 




= 1, 


^(0),[2] 


t 

u 






= 1, 




s 
It 




^(0),[1] 


= 1, 


^(0),[2] 


s 





(2.36) 



A typical partonic cross section requires an amplitude interference, summed over all 
external colors. Such interferences are evaluated in the color bases ( 2.3C| ), ( ^.31| ) as 



riL,L') 



Ai A2A3A4 

where the symmetric matrix CC, 



A2A3A4I A1A2A3A4/ 



A2A3A4 



(2.37) 



c,c'=l 

Ecolo.sTrW*TrM is 
f V -1 N 



OC 



v_ 

N 



-I V N 
N N N"^ 



(2.38) 



One can use table ^ to convert from the A1A2A3A4 helicity configuration label to the label 
h G {1, 2, 3, 4, 5}. The unpolarized partonic cross section is obtained from the helicity sum 



/ > A1A2A3A4 ' 



(2.39) 



A,=±l 



after the usual averaging over initial spins and inclusion of flux factors. For example, 
the helicity sum for the tree-level cross sections for qq gg and qg ^ gq, constructed 
from eqs. ( |2.36| ) and ( |2.28D in either the HV or FDH scheme, i.e., for any 6r, are 



(0,0) 



7(0,0) 

qg^gq 



-2V 



tu 



su 



N 



N 



+ V? 



1 

N 



(2.40) 
(2.41) 



gg were first evaluated through 0{e^) as an interference 



3. One-loop amplitudes 

The one-loop amplitudes for qq 
with the tree amplitude in the CDR scheme |44]. Later they were evaluated as helicity 
amplitudes in the HV and FDH (or DR) schemes [pO| ]. 

Because I^^^ contains terms of order the I^^^\A4^^^)^s. term in the infrared de- 

composition ( p. 11 ) of the two- loop qq gg amplitudes requires the series expansions of the 
one-loop amplitudes through 0{e^). In section 3.1 we present the all-order results in the 
color bases ( 2.30 ), ( ^.31 ), with the normalizations implicit in eq. (12^ 



in terms of integral 
functions whose series expansions are known to the requisite order ||2|, |3| . 

In ref. @ we showed that 0{e) terms in one-loop amplitudes such as M^^q^gg are 
not required for the construction of a numerical NNLO program, once such terms have 
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been subtracted from Ad^g^^gg in the framework of ref. [^]. Thus we need only present 
exphcit formulae for the finite remainders M^g^^gg of the one-loop amplitudes, after ultra- 
violet renormalization ( |2.8| ) and subtraction of infrared divergences ( 2.10| ). We do this in 



section 3.2, for a general 5r scheme. 



In section 3.3 we give a formula for the contribution of the finite remainders A^^^^^^ 



to the NNLO cross section. We compare the result of evaluating the formula in the HV 
scheme to other computations in both the HV and CDR schemes fs^ . 

Actually, the formulas ( [4. 3D and (|4.4| ) for converting the two-loop finite remainders 
■M^^q^gg from one scheme to another are most compactly presented in terms of the 
dependent parts of the one-loop amplitudes at order e; these quantities are presented in 
appendix 

3.1 All orders in e QCD amplitudes 

We now present the renormalized one-loop qq — > gg amplitudes in the color bases ( 2.3C| ), 



( 2.31 ), with the normalizations implicit in eq. (12. Tf ), in a form valid to all orders in e. 

At one loop the crossing properties of the amplitudes are relatively simple, so we 
present the explicit values of the helicity amplitudes for the process qq — > gg. The qg — > gq 
process may be obtained from these by crossing the antiquark into the final state, and 
gluon 4 into the initial state, 

MP'[^](s,t,u) = M(^)'[^](t,s,n), (3.1) 
Mi^^'l'l {s, t, u) = M^^)'''^^ (t, s, u) , (3.2) 
M«'M(.,i,n) = ^M«'['='](t,n,s), (3.3) 



where M^'^'''^^' is defined in eq. (|]2|) with the color bases ( ^l30|) and (|2^ ) using the 



helicity configurations h defined in eqs. ( 2.22| )-( 2.26| ). The factor of s/u in the third relation 



accounts for the permutation of the 5*2 prefactor in eq. ( p. 27 ). The color label c' needed 



for the h = 5 case ( |3.3D is given by c' = 2, 1,3 for c = 1,2,3, respectively. After crossing, 
appropriate analytic continuations are required to bring each function into the physical 
region; for the finite parts for convenience we will give the explicit forms in the different 
analytic regions. 

It is convenient to give the amplitudes in terms of "primitive" amplitudes, which are 
color-stripped building blocks for full amplitudes in any color representation. For the case 
of one-loop amplitudes with a single external fermion pair, the explicit relations between the 



primitive amplitudes and the color decomposed amplitudes were presented in ref. ||3l|]. Here 
we quote the results for the decomposition and then present the primitive amplitudes. The 
same amplitudes, but in a form valid only through 0{e^), as needed in an NLO calculation. 



may be found in refs. 31| 



The first coefficient in the color basis ( p.30|) for qq —)■ gg at one loop is expressed in 
terms of (unrenormalized) primitive amplitudes as [^] 

Ml^)'W(.,i,n) = Afyl^(l+,2,r,3+,4+)-lA^(l+,2r,3+ 4+) + 
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+ iV^A^'[i/2](l+,2„3+,4+), (3.4) 
M«'W(5,i,u) =iV^^(l+, 2r, 3-, 4+) -1^^(1+2^,3-, 4+) + 



+ iV, A^.IV2](i+,2r,3;,4+) - ^Mf 'W , (3.5) 

where, as defined in ref. [^], the L and R superscripts refer to whether the external 
fermion line turns "left" or "right" upon entering a diagram. The "[1/2]" designation on 
the last primitive amplitude represents the subset of contributions with a closed spin 1/2 
fermion. We continue to label the helicities using the "standard" convention with legs 1 
and 2 incoming and legs 3 and 4 outgoing. 

The coefficient of the second color factor for /i = 1 in the basis (2.30) is obtained from 
the first one by permuting kinematics and helicity labels, 



Mf)'[^l(s,t,u) = --Mf)'W(s,u,i), (3.6) 



where the —t/u prefactor accounts for the implicit permutations of the removed Sh pref- 
actors in eq. ( |2.27 ). The coefficient of the second color factor for h = 2 \s given by 



Mi')'['l(5,i,n) =iV^^(l+,2r,3+,4-)- 1^^(1+2^,3+4-) + 



+ iV, A^.IV2](i+,2„3+,4p - ^Mf)'I^] . (3.7) 



For the third color structure in eq. ( |2.30| ), we have 



Mf)'[^l(.,i,n) = ^^(l+,2,r,3+ 4+)+yl^(l+,2r,3+,4+)--^^(l+ 2r,4+ 3+^ 

-^yl^(l+,2r,4+,3+) + 2^A^(l+,4+,2r,3+), (3.8) 

M(')'['l(5,t,n) = ^^(l+,25,3,,4+)+yl^(l+,2r,3-,4+) + ^yl^(l+ 2r,4+,3-) + 

+ ^^^(l+,2r,4+,3-) + 2^^«(l+,4+ 2r,3,-). (3.9) 



where again the ratios of kinematic invariants (including signs) account for the permu- 
tations of the extracted the overall prefactors ( |2.27D from the amplitudes with standard 
ordering of legs (Ig, 2g, 3g, 4g). 

In order to compress the notation a little, we shall suppress the labels of legs in the 
explicit formulas for the primitive amplitudes. We always take them to be ordered (1,2,3,4) 
in the following formulas, and apply permutations as required by eqs. (^)~(|3^). That is, 
define 

aL — aL('\\i 9A2 oAa /iA4\ 

aR — aR (-1X1 9A2 oAs a\4,\ 

^q^lq^2gHgM - ^ l-Lg i-Sg J ; 

aL,[1/2\ ^ .L,[l/2]/-|Ai 9A2 0A3 AXi\ 

q>-lq^2gHg>-i — ^ \^q ^^q ^'^g ^'^g ) ■: 

^q^lg^2qHgH = 2g^ 4^^ ) . (3.10) 
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Then, for example, 2^ , 4+, 3^ ) = g+g-^t^^- 

In this notation, the exphcit values of the independent primitive helicity amplitudes 
— in terms of a set of scalar integral functions — are: 



Ai 



'q^i 9+9 



e(l - eSn) t 
Au 



e{{5-2e)s + 2t) 
(1 - e)(l - 2e)(3- 2e) 



TriW(c 



+ 



e s 



(l-e)(l-2e) 



Tri(^)(i) -sBox(6)(s,i) 



"^q+g-g+g^ 



A 



Q^g 9^9^ 



^q+q-9-9'' 



aR 

^q+q-9-9'' 



A 



^,[1/2] 

q+q~9~9~^ 



A 



g^Q 9^9 



g^q g^g^ 



2(l-e)(l-2e)(3-2e) 
(1 - e6R)st 



q^q g^g^ 
TVi(^)(s), 



■(((3-e)s + t)Tri(^)(s) + 



1 

2^ 



(l-e)(l-26) 
+ ((1 - e)s - t) Tri(4) (t)) - (2 - e)s Box^^) (s, t) 

+ 2i((l - e)s + (1 - 2e)t) Tri^^^tfj + 
+ ((1 - 2e){s^ + 2*2) + (1 - 4e)st) Box(6)(s, t) 



+ 



(1 - e5R)st 



-((l-e)s + (3-2£)t)TriW(t)) 



— es — 2(1 -e)t) Box(6)(s,t) 



+ 



s 

^ 2^ 



(((1 - 2e)s + t) TriW (s) + 2et Tri^ (t)) + 



+ ((l-2e)s-i) Box(6)(s,t) 



= 0, 



e(l - e6R)s 
Au 



{l-e){l-2e) 
+ sBox(^)(s,*) 



(3.11) 

(3.12) 
(3.13) 



(3.14) 



(3.15) 
(3.16) 



((l-e)s + i)TVi(4)(s) + e*Tri(^)(0) + 
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2(1-26) ^^^'^^"^ - t TViW(t) - i(. + 2(1 - 2e)t)Box(6)(.,f) 



(3.17) 



aR 

^q+q-9+9- 



e(l - e5R)s 



((l-e)s + (l-2e)t)Tri(4)(s) 



(l-e)(l-26) 

- et Tri(^) (t)) + (s + 2t) Box^^^ (s, t) 



— ^ Tri(4) (s) - - Box(6) [s, t) , 

2(1 - 2e) ^ ^ 2 ^ ^ ' 



i,[l/2] 



-e(l - e5^j)t 



(l-e)(l-26) 
— TViW(t) + Box('')(s,t) 



(T¥i(^)(s) + Tri(^)(t)) - Box(6)(s,t) 



(3.18) 
(3.19) 

(3.20) 



1 - 2e 



■T>iW(c 



_l_TriW(t) - i((l - 26). + t) Box(6)(., t) . 



(3.21) 



In assigning helicity labels to the above primitive amplitudes we take the quark legs to be 
incoming. 

Here Tri''^-*(s) is the scalar triangle integral in 4 — 2e dimensions with one external 
massive leg, and Box^^^(s,t) is the all-massless scalar box integral in 6 — 2e dimensions. 
The explicit value of the triangle integral is 



TriW(s) = -!I(-s)-i 



(3.22) 



where 



rr 



r(i + 6)r^(i-6) 

r(l-2e) 



with 



n 



(2 



n=l 



TT 



Cs = 1.202057... 



90 



(3.23) 



(3.24) 



In the s-channel (s > 0), the e-expansion of eq. ( p.22| ) is given by using the analytic 
continuation ln(— s) — > Ins — in. The D = 6 — 2€ scalar box integral is completely finite as 
e ^ 0. Its expansions to 0{e^) in the various kinematic channels are given, for example, 
in refs. [0, i. 



3.2 Finite remainders 

Next we tabulate the finite remainders of the one-loop qq — > gg and qg — > gq amplitudes at 
C(e°), defined by M^g^^gg and Ai^J^gg in eq. ( 2.10| ) and color decomposed into Mj^^'^'^'^^^ 
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in eq. ( |2.29D . We write, 



M, 



(l),[c]fin 



h 



-6o(ln(5//i')-m) + %(l-5R) 



(0),[c] 



N 



(3.25) 
(3.26) 



For helicity configuration /i = 1, Bose symmetry under exchange of legs 3 and 4 (t <-> u) 
implies (see eq. (^^)) that 



af^ (s, t, m) 



-— C^l^ (s, M, t) , 

■-6W(s,u,t), 



df^(s,t, n) = (i-l\s,u,t) 



(3.27) 
(3.28) 
(3.29) 



For the h = 1 amplitude, the independent remainder functions a, 6, d, h and j are 



X 1 

6 ~ 4 
1 

'4' 



[1] 



6 



[3] 



0, 



where 



t u 

y = - 

s s 
For h = 2 the functions are 



X = In 



y = In 



an 



1 



(x -y){l- xy) 2 Sx^ - Sxy + lly 



4y3 

(x - y){l - xy) 
2y3 



1 



X 



— 5- ^ X H h 2 + ivr 

X in 

6" + y ' 



X + + 

12y2 4y 2 

6x2 _ ^ ]^-|^y2' 

X 2-y 
2? ^ 



^~yy2_ 5x^ 



4y 

y2 



3y 



3 X 
2y 



+ ivr 



4y^ 



X — y 5 X 



2y 



3y. 



7[2] 



X y 

^2- +171 — 

4y y 2y 

X X 
y + ZTT , 

6y 6y 



(3.30) 

(3.31) 

(3.32) 

(3.33) 
(3.34) 

(3.35) 



(3.36) 
(3.37) 
(3.38) 
(3.39) 

(3.40) 
(3.41) 
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^[3] ^ x{i - 2x) _x^_y;^ 



5 X 5 vr^ 

H X + -Y ^ — 

2 3y 3 2y 



+ 



+ m 



2 — xy + x 



X + -Y 

y 



X Y in 

X \ 

Qy 6 6y 



For h = 3 the functions are 



hi 

.[3 



_5_ 



1 

4' 



_ 1 

6x 
_1 

~4' 
1 

6x ' 
1 _^ 1 

6x Ay ' 
1 

%' 
1 

0, 
0. 



For h = 4 the functions are 



[1] _ (1 - y)(l - ^2 , 6a;2 + ISxy - 2y2 ^ 

*4 ^ 5 ^ I , „ r, A + — 

+ 15xy — 2y'^' 



d 



[2] 



,•[3] 



+ 



(1 -y)(l - xy) 

2y3 



12y2 



[1] _ 



x^ 2 , - y) X , 9 , w 
A H —5 X + \- Z + i-K 

X i-K 

6" + y 



4y 



[2^ 



■X + 



x{2x — y) 

4y2 _ 



,[2] _ 



J2] 
O4 



1 



((X-F)2 + 7r2)- 



4y 

X 

4y 
X + Y 
6y 



X + lOY 3 . 11 
6y 2y 6y ' 



3y 



X(l - 2x) -^2_Y^ XY , TT 



4y2 



+ iiT 



2 y 



^ 5 X 5 



2y Sy 



x{l-2x) ^ x^^ 5 

2y^ y/ y 3yJ ' 



6y 6 6?/ 
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For h = 5 the functions are 



J2] 



J21 
^5 

y[2] 
,[3] 



4y 



2y 6 



— X' + 2 + i^— X 
4y 2y 

X iir 



(1 -y)(l -xy) 



4y 
J__l 
Ay 4 



(X - Yf + TT^ 



12y 



{X-Y)-—Y 



11 

6y ' 



2^ x(l - 2x) 7 1 

" 4y ^^-^) + %-4 



6y 3y ' 

4y V / 2 2y 3y 3 



+ ivr 



5 



6y 6 6y 



(3.60) 
(3.61) 
(3.62) 

(3.63) 
(3.64) 
(3.65) 

(3.66) 
(3.67) 



For the HV scheme {5r = 1), the results ( |3.25| )-( |3.67| ) for the finite remainders of the 
one- loop helicity amplitudes are in complete agreement with those of ref. [42|. 

3.3 Comparison with CDR results 

Results in the CDR scheme are usually phrased in terms of amplitude interferences, summed 
over all polarizations and colors. In the NNLO cross section, the one-loop amplitude en- 
ters interfered with itself. The contribution of the one-loop finite remainders to the NNLO 
qq — > gg or qg gq cross section, summed over all helicities and colors, is given by 



j(l,l)fin ^ 

A,=±l 



V^' Ai A2A3A4 l-'^'^Ai A2A3A4/ 



(3.68) 



Using the color sum matrix OCij in eq. (2.38), and table |l]to relate A1A2A3A4 to h helicity 
configurations, the color and helicity sum in I^^'^)^'^ may be evaluated in terms of the above 



explicit expressions (p5|)-(p7|) for M^^^'''^^' , in a general 6ji scheme 



In ref. |g] , the one- loop-squared contributions to the qq — > gg and qg — > qg cross sections 
were given in the CDR scheme, but not in a convenient form for comparing to our results, 
because the infrared poles were not organized in the same way. (Also, the terms related to 
ultraviolet renormalization in eq. (4.1) of the original version of ref. H] require correction.) 
However, the authors of ref. |^] have kindly supplied us with their versions of the values of 
eq. ( 3.68| ) in both the HV and CDR schemes. These expressions agree precisely with our 
result for eq. ( |3.68| ) in the HV scheme {6r = 1). 
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4. Two-loop QCD amplitudes and finite remainders 



We generated the Feynman graphs for qq — > gg using QGRAF |54], from which a MAPLE pro- 
gram was constructed to evaluate each graph. As a cross-check, some of the diagrams were 
evaluated using FORM [^]. We employed the general integral reduction algorithms devel- 
oped for the all-massless four-point topologies [0, 18, 16, in order to reduce the loop 
integrals to a minimal basis of master integrals. To put the integrands into a form suitable 
for applying the general reduction algorithms, spinor strings were converted to traces over 
7 matrices, by multiplying and dividing by appropriate spinor inner products constructed 
from the external momenta. Evaluating the traces then gave dot products of momenta; 
any terms containing an odd number of Levi-Civita tensors ef^'^'^P vanished upon integra- 
tion. The gluon polarization vectors of definite helicity can be incorporated using some 
minor extensions of the integral reduction techniques |^]. Here we incorporated the gluon 
polarization vectors in a slightly differently fashion than described there. When forming 
the traces over 7 matrices we included also the polarization vectors expressed in terms of 



spinors |32], following the methods used in, for example, refs. |35, 56|. Because the polar- 
ization vectors are four-dimensional objects, they distinguish between 4-dimensional and 
(— 2e)-dimensional components of the loop momentum. We evaluated the resulting inte- 
grals, containing (— 2e)-dimensional components of the loop momentum, using the methods 
of ref. p. (See refs. |^, |ll|, |4^] for alternative multi-loop helicity techniques.) 

After all the tensor loop integrals in the amplitudes have been reduced to a linear 
combination of master integrals, the next step is to expand the master integrals in a 
Laurent series in e, beginning at order 1/e^, using results from refs. [^, |l^, [l^, 18|. It 



is straightforward |57| to express the results solely in terms of polylogarithms 



Unix) 
Li2(a:) 



EX 

-It-" 



— Li„_i(t) , 



t 



t) 



(4.1) 
(4.2) 



with n = 2, 3,4. The analytic properties of the non-planar double box integrals appearing 
in the amplitudes are somewhat intricate ||2|, Q ; there is no Euclidean region in any of the 
three kinematic channels, s, t or u. So we do not attempt to give a crossing-symmetric 
representation, but instead quote all our results in the physical s-channel (s > 0; t, u < 0) 
for both the qq — > gg and qg gq kinematics, eqs. (|2.l| ) and (^.21) . 



4.1 Finite remainders 



The two- loop finite remainders are defined in eq. ( 2.11| ) and are color decomposed into 
j^^^2),[c]fin ( |2.29| ). Their dependence on the renormahzation scale color factors N 

and A^, and scheme label may be extracted as 



M, 



(2),[c]fin 



+ (2Rq + 2Rg + bQQ^^^\\n{s/n 



zvr -I- 



-^-CA + \cF)CF{l-5Rf + 



1 



m) + -CFQQiTi]{l-5B) 



^(0),[c] ^ 
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+ 



(l),[c]fin 
h 



+ 



M, 



CaV. 



(1 - <5ij) + 



Ar2 



(2),[c]fin 



-26o (ln(s//i2 



ivr) +qS^^)(1 



+ (5o 



h 



CaV, 



+ — /' 



(1 



c = l,2, 

^(l),Wfin ^ 



(4.3) 



(4.4) 



The ^-dependence is a consequence of renormalization group invariance. 



The tree and one-loop functions, M, 



(o),W 



and M, 



(l),[c]fin 



h 



are given in eq. ( |2.36| ) and 



eqs. ( 3.25 )^( 3.67 ), respectively, while 6o and bi are given in eq. ( |2.6| ). The following com- 
binations of color constants also appear in eqs. (E^) and ([4.4[), 



1, 



+ -TrNj 



(99) 



5c. 

6 

D 2 



Rn 



V192 
5 



864 



VT" 



+ 



16 



(4.5) 
(4.6) 
(4.7) 
(4.8) 



where the combination AR„ appears as well in the 6r dependence of the gg — > gg amplitude. 



eq. (5.10) of ref. |]. The quantities 



(l),[c]e,5fl 



are the (5R-dependent parts of the 0{e^] 



coefficients of the one-loop amplitude remainders, after subtracting the poles according to 
eq. ( p.lOj ). The explicit values of M^^'^^'^'''^^ are tabulated in appendix ^. Finally, the 



quantities 
h = I and h 
amplitudes ^ 



are only non-vanishing in the case of the simpler helicity configurations, 
3. They seem to be related to the 5r dependence of one-loop splitting 
Their explicit values are given by 



(1),[2] 



X 

6 



(1),[3] 



0, 



V, 



(1),[3] 



0. 



(4.9) 
(4.10) 



X 

6 ' 

6x ' 6x 
The coefficient functions A, B, C, D, E, F, G, H, I, J, K, L depend only on the Mandelstam 
variables. In appendix we give the explicit forms for the independent finite remainder 
functions appearing in eqs. ( |4.3| ) and ( [4.4p . 

We have compared our results for the independent two-loop finite remainder functions 
^^2),[c]fin ^-^j^ corresponding results obtained contemporaneously in ref. [42|. The results 



agree completely, once a correction is made for a slightly different definition of H^'^\e) in 
eq. (|2.18| ). (We always dress the 1/e pole with (;U^/(— s))^*^; ref. [^] sometimes dresses it 
with (^7(-t))2^ or (^2/(-^x))2^) 
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4.2 Comparison with CDR results 



Finally we discuss conversion from the HV scheme results reported in section Q to the 
CDR scheme used in ref. In the CDR scheme, one usually computes the interference 
of amplitudes, summed over all external colors and (2 — 2e) polarizations. The generic 
one-loop/tree interference encountered at NLO is 



2 Re I^f^ = 2 Re 



color, hel 



(4.11) 



Inserting the infrared decomposition ( |2.10| ) for Ain^ into eq. ( [4.11 ) gives 



rim 



2Re 



(0)\ 



?(l,0)fin 

"T -'rs 



color, hel. 



(4.12) 



where 



(o)\ 



color, hel. 



(4.13) 



It is well-established from explicit calculations and general arguments |46, |50|, 51 1 that 
the finite remainder ( 4.13| ) has the same value in the HV and CDR schemes, in the limit 
e — > 0. Essentially, the treatment of unobserved partons is the same in both schemes, 
so the infrared divergences should take the same form, when expressed in terms of the 
lower-order-in-CKs amplitudes. 

It is natural to expect the same pattern to hold at two loops. (Indeed it does for the 
99 ~^ 99 amplitude The two- loop/tree interference is 



2Re/S°) 



2 Re [i^^n^\^ 



(0)\ 
n I 



color, hcl. 



2Re Y + 



color, hel. 



(4.14) 



+ (4.15) 



where 



I, 



(2,0)fin 



2 Re Y [{M^n^^^\M^:, 



(0)\ 



(4.16) 



color, hcl. 



Note that /^^^ and J^^) are the same operators in the HV scheme as in the CDR scheme. 

We have interfered the color-decomposed finite remainders of the two-loop qq — > gg 
and qg — > gq helicity amplitudes in the HV scheme, as given in section 4.1, with the tree 
amplitudes given in eq. ( 2.36[ ), summing over all external helicities and colors with the help 
of eq. ( 2.3^ ). This sum gives precisely the same result as the corresponding quantity ( 4.16| ) 
in the CDR scheme, as evaluated in ref. ||5|, after accounting for the slightly different 
definition of iJ^^^ that we used in eq. ( ^.18 ). This result provides additional evidence 
that eq. ( [4.16 ) should be the same in the HV or CDR schemes for general two-loop QCD 
scattering amplitudes. A similar conclusion has been reached independently [p2|. 
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5. Amplitudes in pure = 1 super- Yang-Mills theory 



QCD amplitudes may be converted easily to amplitudes in = 1 pure super- Yang-Mills 
theory by modifying the fermions to be in the adjoint representation and by altering their 
multiplicity to be a single Majorana fermion. This theory is then supersymmetric with the 
fermion the gluino superpartner of the gluon. Besides the inherent interest in supersym- 
metric theories, a practical consequence is that supersymmetry imposes a set of powerful 
identities that provide non-trivial checks on the amplitudes, including their finite parts. 
The supersymmetry identities have been applied previously to the same one-loop ampli- 
tudes discussed here [|^], but only through 0{e^), as needed in an NLO calculation. Here 
we extend the one-loop discussion to include all orders in the dimensional regularization 
parameter e, which are relevant at two loops via Catani's formula ( |2.11| ). Then we verify 
the identities at two loops. 

5.1 Supersymmetry Ward Identities 

One set of identities implies that "maximal helicity violating" amplitudes vanish for any 
supersymmetric theory and any number of loops, 

M^''^''{gt,9^,gt,...,g+) = 0, (5.1) 
M^''''^igt,g2,9t,---,9t) = 0, (5.2) 

where g and g denote a gluon and gluino, and the superscripts denote helicities. In this 
paper we use the convention that legs 1 and 2 are incoming and the remaining ones are 
outgoing. Other identities relate the non- vanishing supersymmetric helicity amplitudes for 
external gluons alone, to amplitudes where some of the gluons are replaced by gluinos. In 
particular, the two- gluino two- gluon amplitude can be expressed in terms of the four- gluon 
amplitude: 

M^''^^i9t,g2,9s,gt) = ^M^''^^{gt,g,,gs,gn (5.3) 



This identity is somewhat more stringent than the ones in eqs. (5T) and ( ^.21) because it 
relates distinct amplitudes containing the most intricate infrared divergences, up to order 

poles at two loops. These relations are crossing symmetric, when a crossing symmetric 
definition |Q of the spinor products is used. 

For the four-point case, the pure-gluon identities (5J) have already been checked at 
one loop to all orders in e, and at two loops through 0{e'^) in the FDH regularization 
scheme [|4|,||. 

5.2 Color and infrared structure 

Since the gluinos are in the adjoint representation we use the same color basis as used for 
the four-gluon helicity amplitudes 

= ShxY, Tr'^] X Mf , (5.4) 

c=l 
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where 



Tr 



tr (T'*^ 2^«2 2^«4 "J 
tr (T'*^ j^f2 2^04 2^13 ^ 
tr (T'*^ j^a4 2^(12 j^fia ~j 
tr (T'*^ j^aa 2^<i2 jifi4 "J 
tr (T'*^ j^«3 JT14 ji«2 ^ 
tr (T"^ j^f 4 jTia 2^12 ~j 
tr(r"iT"2)tr(r'^3 2-a4)^ 

tr(r''ir"3)tr(r"2r"*) , 
tr(r''iT"*)tr(r''2r"3) . 



(5.5) 



A reflection identity implies that the c = 4, 5, 6 coefficients are equal to the c = 3,2,1 
coefficients (respectively), so there are really only six different coefficients for each h, namely 
j^{L),[c]^ c= 1,2,3,7,8,9. 

Using C, P, and T, for the case of two external gluons and two external gluinos, there 
are five independent processes to consider, including processes related by crossing: 



h = 


1 


g{pi 


+)+9{P2,-) - 


g{p3, 


+) 


+ 9{Pi 


+)> 


h = 


2 


g{pi 


+) + g{p2,-) - 


g{p3, 




+ g{p4 


+)> 


h = 


3 


g{pi 


+) + gip2,-) - 


g{p3, 


+) 


+ g{p4 


+)> 


h = 


4 


9{Pi 


+) + giP2,-) - 


g{p3, 




+ 9{P4 


+)• 


h = 


5 


g{pi 


+)+giP2,+) - 


g(.P3, 


+) 


+ g{p4 


+) • 



(5.6) 
(5.7) 
(5.8) 
(5.9) 
(5.10) 



The latter three processes are obtained from the first two by crossing. However, just as for 
the quark case we keep them distinct, because the crossing properties at two loops are in 
principle nontrivial. 

The infrared divergence structure is similar to that of gluon-gluon scattering ampli- 
tudes 0| . For the case of = 1 pure super- Yang- Mills theory, in the basis (^]B) the matrix 
!<•> isl, HI 



£V(1) 



r(i~e) W 



Ne 



(5.11) 









(S-U) 


V (T-U) 



(T-U) (S-U)N^ 

7V(S + U) (U-T) (s-T) 

iV(T + U) (T-S)(U-S) 

iV(T + U) (T-S)(U-S) 

iV(S + U) (U-T) (S-T) 

iV(S + T) (T-U) (S-U) 

(S-T) (S-T) (S-U) 27VS 

(U-T) (U-S) (U-S) (U-T) 2m 

(T-S) (T-S) (T-U) 2NT ) 



-22- 



where S, T and U are defined in eq. (l2^ ). For = 1 super- Yang- Mills theory the first two 
coefficients of the /3-function are 



bo = -jCA , hi 



The -T*-^-* operator for super- Yang-Mills theory is 

f{2) 



3^2 



(5.12) 



'FDH 



I (e,At;M) + — 



+ 



+ 



g+ev(i)r(i _ 2e) / 6o 



rfi 



— -t- -n-FDH 



2e,;U;M) + 



+ -H"FDH(e,/^;M) , 



where 



7V-SYM 

-"■FDH 



-H'FDH(e>^;M) 



4er(l 



and 



SYM 

g ;FDH 



2e 



2 



4(/f(2)) 



SYM 1 I tJ- 

g /FDH + 



(2) 



(5.13) 

(5.14) 
(5.15) 



vr 
16 



(5.16) 



The equality of (-ffg^OpDH (-^^ Ofdh ^ consequence of supersymmetry. Equa- 



'g Jfdh 

r(2)^SYM „„J /'fr(2)\SYM 
9 JfDH "^^^^ \^g JfDH 

tions ( |5.14D ~( ^T6|) are obtained from the QCD formulas, eqs. ( |2.17| )~( ^3o| ), by the re- 
placements 5r — > 0, Cf — > Ca and TrNf CaI'^ for converting to a single adjoint 
fermion in the FDH scheme. The operator H. defined in eq. (2.21) does not explicitly 
depend on the fermion representation. 

(0),[c] 



The tree amplitudes in this color basis are given by , where 



M\ 



(0),[c 



(0),[c 
h 

(0),[c 



(0),[c 



,>(0),[c] 



0, 



U 



c= 1,2,...,9, 

c = 7,8,9, /i = 1,2,3,4,5, 
c=l,6, /i = 2,4,5, 
c = 2,4, /i = 2,4,5, 
c = 3,4, /i = 2,4,5. 



(5.17) 



5.3 One-loop amplitudes in pure = 1 super- Yang-Mills theory 

We now present the results for one-loop two-gluino two-gluon scattering in a format valid 
to all orders in e. By comparing these results to the corresponding = 1 pure Yang-Mills 



four-gluon amplitudes ||9|] , we then verify that the identities (|5.2|) and (l5.3|) do indeed hold 
at one loop to all orders in e when the FDH scheme is used. 

To check the supersymmetry Ward identities at one loop, we use a crossing symmetric 
representation of the amplitudes directly in terms of one-loop scalar integral functions. 
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Thus it is sufficient to explicitly present only the h = 1,2 helicity cases. Using formulae 
from ref. [^] for obtaining the = 1 supersymmetric amplitudes from the primitive 
amplitudes, the coefficients of the first color structure Tr^^l in eq. ( ^ ) for h = 1,2 are 
given by 

M['^'^'\s,t,u) = iV(^^i(i+,2r,3+,4+)+A«(l+,2r,3+,4+) + 

+ ^^'[1/21(1+ 2r, 3+4+)), (5.18) 
Mi^)'W(5,t,u) = iv(A^(l+,2^r,3;,4+) + .4^(l+,2r,3-,4+) + 

+ A^'[V^l(1^2„3„4+))-^Mf'W. (5.19) 



e 



For the second color structure in eq. (p.5|), we have 



Mi'^'^'\s,t,u) = -ivl(^^(i+ 2gr,4+ 3+) + ^^(l+,2,-,4+,3+) + 



+ ^^'[i/2l(l+,2r,4+,3+)), (5.20) 
M«'['l(s,t,n) = iV^(^^(l+ 2r,4+,3-)+yl^(l+,2r,4+ 3,-) + 

+ ^L,[i/2](,+^2r,4+,3p) - ^Mf . (5.21) 

The third color configuration Tr''^' has coefficients given by, 

Mj^^'f^^(s,t,u) = 0, (5.22) 

Mi')'['](s,t,n) = 2Ar-A^(l+ 37,23,4+) - ^Mf^'P' . (5.23) 

In the above, hQ is the first /3-function coefficient for = 1 super- Yang-Mills theory. The 



remaining color coefficients are given in terms of these |61, 35 1, 

M(^)'W(«,t,n)=M«'[^l(.,t,n), 
Mi^)'['l(s,t,n) =M«'['l(s,t,n), 
M«'[^](s,t,n) = M«'W(.,i,n), 

m(^)'['] (s, t,u) = ^ [mI^^'^'^ {s, t, u) + m(^)'['1 {s, t, u) + iV^ W is, t, u)] 

M^^^'f^^ (s, t, u) = M^^^'f^l (s, t, -u) = M^^^'f^l (s, t, u) , (5.24) 

for all helicity configurations h. 

Inserting the explicit values of the primitive amplitudes and taking the FDH scheme 
{6r = 0), yields for the three independent color factors of the h = 1 helicity configuration 
in eq. (|5^), 

(5.25) 
(5.26) 
(5.27) 



M«'W(. 


t, 


u) 


= 


<)'P1(. 


t, 


u) 


= 


<)'[^1(. 


t, 


u) 


= 
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For the h = 2 helicity configuration (p. 



{s,t,u) = N 



s 2n + e(4s + t) 



+ 



2u 
1 

2u L 



TriW(.)-±'" + ^('^ + ^^ TViW(t) + 
2u 1 — 2e 



1 - 2e 

2^2 + + 2*2 - e(4s2 + 5ts + At^)] Box^^) (s, t) 



,>(1)'[2] 



M 



e st 



1 - 2e 2m 



TriW(s) _tTVi(^)(u) + 



+ 



,>(1),[3] 



M 



2u 
e st 



s,u) 



^^(0),[2] 



(5.28) 



(5.29) 



+ 



1 - 2e 2n 
s(2s + e(4it + t)) 
2u 



TriW(t) -sTri(^)(n) + 



Box(6)(t, 



^^(0),[3] 



(5.30) 



For helicity h = 1, the amphtudes clearly satisfy the super symmetry Ward iden- 
tity ( |5.2| ) exactly. To check identity (^I^), we compared the amplitudes in eqs. ( 5.2^ ), 
(5.29) and ( p. 30 ) to the corresponding four-gluon amplitudes. The latter amplitudes may 
be obtained from section 3.1 of ref. by letting Nf N and 5o ~^ ^0 in eq. (3.1) of that 
paper. (Note that the 'W = 1" there refers to a chiral multiplet consisting a scalar and 
a fermion, while here we are considering a pure super- Yang- Mills multiplet consisting of a 
gluon and a gluino.) The result of this comparison is that the supersymmetry identity ( |5.3| ) 
holds to all orders in e when the FDH scheme is used. 

We note that the symmetry relation evident between eqs. ( 5.2g| ) and ( 5.30| ), 

M^^^''^^ {s, t, u) = M2^^^''^l {t, s, u) , (5.31) 

which involves a swap of a gluon and gluino leg, is a consequence of the Bose symmetry 
of the corresponding four-gluon amplitude, because the gluon legs being exchanged in the 
latter amplitude have the same helicity (when both are considered as outgoing states). 

After performing the subtraction of infrared divergences using eqs. and dsll] ), 

we obtain 



(l),SYM,[c]fin 



h 



-6o(ln(.//x2) - .vr) + ^-n] + AT a^'^^^ , 



(5.32) 

c = 1,2,3, 

- yh ^ c = 7,8,9, (5.33) 

where the one-loop supersymmetric remainder functions are given in terms of the QCD 



M, 



(l),SYM,[c]fin 



SymJcI 
9h 



ones. 



symJcI 

symJcI 
9h 



[c] _ f\c] 

h "h 

SYM,[1] 



SYM.[2l 



SYM,[3l 



c= 1,2,3, 
c = 7,8,9. 



(5.34) 
(5.35) 



The A^M^'^^''^^/2 term in eq. ( |5.32| ) is a consequence of the finite shift between the FDH 
scheme used in super- Yang-Mills theory and the HV scheme used in QCD. 
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5.4 Two- loop amplitudes in pure = 1 super- Yang-Mills theory 

An immediate consequence of the one-loop supersymmetry identities holding to all orders 
in e is that all infrared divergent terms at two loops also satisfy the same identities. Since 
neither I^^^ in eq. ( ^.11 ) nor I^'^^ in eq. ( ^.13 ) depend on whether the external lines are 



gluons or gluinos, the Catani formula (|2.11| ) dictates that divergent terms at two loops 
must satisfy the same supersymmetry Ward identities as the tree and one-loop amplitudes 
do. 

This then leaves the question of whether the two-loop finite terms satisfy supersym- 
metry Ward identities. For the case of four-gluon scattering, the finite remainder terms, 
after subtracting the Catani terms (including their finite parts), have been presented pre- 
viously p. We have carried out the analogous subtraction using eqs. ( p. 11 ), ( 5.11 ) and 



(|5.13D . The result of this is, 

(bof (ln(s//i2) - ivr)2 + 6i (Hs/fi^) - in) 
26o (ln(s//x2) - ivr) m^^)'^™'!^'^" + ^^'^ + ^'t^l 



^(2),SYM,[c]fin 
h 



^(2),sYM,[c]fin ^ (ln(s/;"') - iTT) M^^)''™'!^''^'^ + iVG^""'[^] 



c = l,2,3, (5.36) 



c = 7,8,9. (5.37) 



where we have verified that the finite remainder functions match those of the pure gluon 
case, i.e. 

^SYM,[C] ^ ^SYM,[C] ^ ^5_3g^ 
^SYM,[c] _ ^SYM,[c] _ ^SYM,[c] 

^SYM,[C] ^ _i;^SY_M,W (5_40) 

where X G {A, B, G}. The functions x^_^j}-^ and x^_^_^^ for gg — > gg scattering in pure 
= 1 super- Yang- Mills theory are given in ref. Q. (Note that in that reference an all 
outgoing definition of helicity is used.) The relations (|5.38|) - (|5.40|) are precisely equivalent 



to the content of the supersymmetry Ward identities |33], after removing overall factors 
and the divergent terms, and extracting the N and n dependence. Thus they provide a 
direct, nontrivial check on the finite remainders. 

Because the adjoint color indices of the gluino fields are identical to those of gluons, 
and only structure constants /"^'^ appear in the two-loop Feynman diagrams, the color coef- 
ficients for two-gluino two-gluon scattering obey the same group theory relations identified 
for gg gg in ref. 0, 

^SYM,[9] ^ 2(^f "'[^1 + + - i^f "'[2] , (5.43) 
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and 

^SYM,[3] ^ _^SYM,[1] _ ^SYM,[2] ^^-^ 

6. Conclusions 

In this paper we have presented the two-loop ampUtudes for quark-gluon scattering in 
QCD and gluino-gluon scattering in = 1 super- Yang-Mills theory, including the full 
dependence on external colors and helicities. We confirmed that, as in the case of gluon- 

^ (2) 

gluon scattering, there is an additional 1/e pole term, H in eq. ( p.21D having nontrivial 
color dependence. This additional term vanishes after interfering it with the tree amplitude 
and summing over colors. We investigated the dependence of the amplitudes on the variety 
of dimensional regularization employed. The FDH scheme respects the supersymmetry 
Ward identities, in what represents the most detailed test to date of these identities, and 
of the FDH scheme. The scheme dependence of the quark-gluon amplitudes is much more 
intricate than was the case for the four-gluon helicity amplitudes [^. 

The two-loop QCD results, when interfered with the tree amplitude, summed over all 
external colors and helicities, and converted to the CDR scheme, are in complete agreement 
with the previous results of Anastasiou, Glover, Oleari, and Tejeda-Yeomans [^. We also 
expressed the one-loop-squared contribution to the NNLO qq gg and qg — > gq cross 
sections in terms of one-loop finite remainders. These also agree with suitably subtracted 
expressions m the CDR scheme p|. 



So far the new two-loop amplitudes have been implemented in only a handful of phe- 
nomenological studies We may anticipate that once general algorithms for dealing 

with infrared divergent phase space integrations at next-to-next-leading-order are devel- 
oped (see e.g. ref. [^^), many more phenomenological applications will follow. These 
applications would include the implementation of the two-loop amplitudes of this paper, 
or those of ref. Q, as ingredients in a numerical program for computing dijet production 
cross sections at hadron colliders at NNLO in QCD. When this task is accomplished, the 
intrinsic precision on the QCD predictions should reach the few percent level, providing a 
stringent test of the Standard Model at short distances. 
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A. Finite remainder functions for QCD 

In this appendix, we present the explicit forms for the independent finite remainder func- 
tions for the processes qq — > gg and qg — > gq in QCD, which appear in eqs. (|4.3| ) and (0). 
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For the helicity h = 1 configuration in eq. ( 2.22 ) and color factor Tr'^^ in eq. ( 2.3C| ), 
the finite remainder functions are: 



[1] 



[1] 



D 



[1] 



[1] 



1 



48y2 

1 2'5 47'5 

-(180x2 + 239xy + 99y'^)X + —x + 



lUy 



864 288 



y 



ITT 



1 



24y' 



+ 



■(30x^ + 51x^2/ + 22xy2 _ 3y3)j^ _^ 



144y 



(180x2 + 239xy + 99^^ 



1 2 r ^ 3, 115 139 

(4x^ + 5xy + 7y^)X (x - y)Y H x H w + 

g V y y 7 8^ ^ 72 72 ^ 



+ in 



[1] _ y^ - 2x- 



8y2 8x 2y 



16y2 



16x 



16 



(X-y)2+7r2 + 



+ T^(4x2 + 13xy - 3y2)X - ^(x - y)y - II + 
Iby 8 32 



+ ivr 



y^ — 2x^^^ x — y^, 3 — x 



8x 



16y 



37 



24y2 
+ m 



(5 - xy + x")^" + —(9 + 5xy + 27x2)X + — x - ^ 



12y 



r(5 - xy + x^)X + - — (9 + 5xy + 27x 



72y- 



[1] _ x(2x-i)^2 , i-2y 



8y2 



"{x-y) / 



8. + + 

19 



1 3 7 

+ -^{y^ + 7xy - 6x2)X - -(x - y)Y + — x - — y + 



+ ITT 



x(l — 2x) 1 — 2y^^ 6x — 7y 
^ -X H —Y + ^ 



4y2 



4x 



8y 



X 5 ivr 

— A X H X , 

36 54 36 



=0. 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



(A.5) 

(A.6) 
(A.7) 



Bose symmetry under exchanges of legs 3 and 4 (t ^ u) implies that for h = 1 the finite 
remainders of the second color configuration in eq. (|2.30D can be expressed in terms of the 
first ones, i.e. 



Z^^\s,t,u) = --ZPis,u,t), 



(A.8) 



where Z £ {A, B, C, D, E, F, G}. 
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For /i = 1 in eq. ( |2.23| ) and color factor Tr^'^^ in eq. ( 2.30| ) 



H 



[3] 



r[3] 



[3] 



[3] 



L 



[3] 



48y2^ 18y^ 



x-y , 
8y 



{l-xy)({X-Yf + TT^) + 



+ ITT 



24y 



r(-6-25x + 7x2 -4x^)X 



24xy 



-6 - 25y + - 4y^)Y + 



7(x - j/) 
6y 



1 — X^^n 1 — y 



16y 



16y 



x-y , 
8y 



(l + 2xy) + ^ 



+ ^(l + 12x)y-^(l + 12y)X + 



1 — X, 



-Y 



x-y 



8y 8y I2y 

x(2x2 + X + 5) 2 V + y + 5 2 , ^ .oq^ 

2V + — (39x- 55y)X- 

- ^(392/ - 55x)y - ^x(x - y) ((X - y)^ + + 



r(2x2 + X + 5)X 



2y^ + y + 5 



12y2 

x(l-2x)^2 l-2y,,2 1 



Y 



13(x - y) 



8y2 



-X' 



8x 



12x " 24y 
xix-y)({X -Yf +Tr^) + 



+ ^(17x - 19y)y - ^(17y - 19x)X + 

'x(l-2x) 1-22/^ 19(x-y) 
zvr ^ — A ^ Y - 



4s' 



4x 



24y 



(A.9) 



(A.IO) 



(A.ll) 



(A.12) 
(A.13) 



For /i = 2 in eq. ( |2.23| ) and color factor Tr^^^ in eq. (|2.30| ) 



A 



[1] 



2y3 



(x3Li3(-x) - Ux^ + 2/)) X - -(3x3 ^ y3)Li4(-x) ■ 
- ^(x - y){l - xy)X^U2{-x) - j^{7x' - ?,y^)X^ + 



2 4 

+ y (2x3 - y^)X^ + Y^(24x3 + 952/^) 
^7r2Li2(-x) + ^Xy[x^ - 27r2) + 



+ 



144y3 



(-24x2y - 77y3 + 21xy2 + mx^)X'^ + 
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+ (576y + 401y^ + 68V + 172)7t^X 



+ 



+ 



48y2 



(352/2 + _ Sxy){2U2{-x)X + X^Y - 2Li3(-x)) + 

1 



+ (12xy + 241y2 _ 1\x^)Q^ 



72y 



:(452/^-49xy + 86a:^)7r^- 



288y3 
1 

86V 



(2850x2 + 609xy - 2273?/2)X 



475x 36077 



288y 3456 



+ 



Li.(Li3(-)-C3)-^X3 + 

+ Y^(y - ^)(l - xi/)X(3Li2(-x) - 7r2) + + 

+ Jyx2-C3-^(5x2 + 20x + ll) + 

+ 7^(-24^^y - 7V + 21xy2 + 86x=^)X2 + 
48^* 

+ 7r?^(35y2 ^ g^2 _ ^xy){U2{-x) + Xy) - 



24y2 
1 



144y 



(475x3 + 99x2y + 279xy2 - 505y3)X + 



+ 



864y' 



r(32x2 + 5155x + 2273) 



(A.14) 



B. 



[1] 



1 

VL 2 



(4 - 5x - 5x2)Li4(-x) + (5x3 ^ ^.2 _ 4^ ^ l)XLi3(-x) 



TT 



(llxy + x^ + 2y^)C3X + y (17xy2 + Igx^y - + 2x3)Li2(-x) + 
1 vr^ 

+ — (9y + lly2 + 2y3 + Q)X^ + _(3a;3 + lla;y2 + 4y3)^2 _ 



TT 



1440 



(112x3 + 1128x2y + 916x/ + 45y3) 



+ 



y 



^(1 + 13x)Li4 (^-^ + ^(12x2 + 15xy + 7y2)Li4(-y) - 

- (1 - 5x)(^XLi3(-2/) - ^y(Li3(-x) - Cs)) + 

+ ^(3x2 + 2y2 + 3xy)X2Li2(-x) - J(31xy + 9y2 + 24x2)x2y2 _ 
8 8 

- 2y2rLi3(-y) - (2 - x)XyLi2(-x) - -^(-3y2 + 12x2 + llxy)y^ + 



48 



1 



+ xyLi2(-x)F2 + _(6;r2 + llxy + 10y2)yx3 - 
- ^(2 - j/)(3x - y)Xy3 _ l(2ar2 - 2x - I)7r2y2 + 



-30- 



1 

4y3 _ 



2y 

(4a;y + Sx^ + 9y^) ( Li3(-x) - XLi2(-x) ) + 



+ :^(720xy2 - lly^ + SQGx^y + 108x=*)C3 - 
6b 

- -^(1202/2 + 60y3 + 43 + TSy)^^ - 
18 

7^2 -I 

- — (gix^ + UOxy + 160y^ - 132)X + 

18 

+ i (^3Li3(-y) + 3Li2(-x)y - ^{lOx + 3)XY^^ - _ 
^ (416x^ + 27x^ + 742^2 + 190 + 902x)X'^ + 



144y 



2 



+ ^(191 + 27x)y2 - 2^(463 + 366x - 35x'^ - 54x^) + 



^(54x3 + 17a;2y + lOOy^ + 3Qxy^)XY - ^^^^^ 



^ (291y2 + 305x2 + 383xy)X-^ + Hl^ + 



72y2^ ' ■ 72y ' 2592 

+ ITT 



\-y{7xy + 3y^ + 6x^)Li3(-y) - ^(6^ + x)Li3(-x) 

1 TT^ 

- -(1 - 3y - 3y2)C3 - — (2y2 + 8xy + x^)X - 

- ^ (7(3 + - y)xu2{-x) + 

y^ \A 

+ {3xy + 3x2 ^ 2y2)Li2(-x)y + 

1 TT^ \ 

+ -(6x2 ^ ^ 8^y)xF2 _ _(4^2 + + 9)y j + 
+ ^ (^3y - 4)YX^ + 2Xr2 - ^(x - y)Y^ + 2Tr'^X^ + 
+ ^ (^-Li2(-x)(3x2 - 3y2 + 4xy) - 

- ^(34 + 42y + 33y2 + 40y^)X^ + 

+ -^(83x + 47x2 + x^ + 49)7r2^ + 
+ ^(29y' + 30x2 ^ ig^y^xY + ^(11^ - ^y)Y^ + -^7r2 - 



- 31 - 



72y3 
1 



(-504?/ - 19V - 14%^ - 359)X + 



+ ^7r^(-91y + 5Ax)Y - ^-^(48?/^ + 40?/ + 61) 



72y 



^[(2.2-8x + 2)(Li4(-?/)-Li4 (-^)) + 

+ (2x2 -8x- 5)Li4(-x) - ^{2x^ - 8x - 3)Tr^X^ + 



TT 



+ (1 + 4x - x') ( 2(Li3(-a;) - Cs)^ - Y^'^i-"") ) " 



1. 



1 



12 



(1 - 4a; + x')Y\^—Y - 2(Li3(-x) - Cs) - ^Y'X + -Y 



+ -^(lla:;^ + 12a;2 + 117x + 24)7r^ 
180 



+ 



+ 



1 

8?L' 



3(2x2 + y2 - 2x?/)XLi2(-x) + 12(1 - x'^)U3{-y) 



- 3(5x2 + l + 4x)Li3(-x 
+ (4y2 + _j_ 43.^)7^2 



(29x + 41x2 + 15x3 + 9)C3 + 



+ 



+ 



2y2[ 



3(x - l)yLi2(-x) - ^(1 - x)Xy2 + ^{y- 2)X^Y ■ 



^(7x-3)7r2F + l(9-13x)X=^ 



+ 



+ T^TT^i-'^^x^ - 6x^ + 5 - 18x - 22x2)^2 - 
32y^ 

3_ ( 12x(3x2 + 5x + 8)XY + (15 - 118x - 59x2 _ i8a.3N^2\ 

96y^ \ ) 

- -^(2x + 3x2 ^ y^^y2 ^ J_(i2x2 + 19x - 39)X - 
16?/ 32?/^ 

3 , ,,,, 187 17x 

(u- 1)1^ H h 

8?/^*^ ^ 128 32y 



+ m 



^ f8Li3(-x) - 2XLi2(-x) - X3 - 8C3) + - - 

13 1 
— (x2 + 4x + 9)^2 + _(a;2 + 4x + 5)Li2(-x) - -^-k^ 



16y 
1 



+ T7^(4x2 - 2x + 5)X - -^(5y - 6)y 



1 



16y3^ ^ 8y2^ ' 32y2 



(17x + 63) 



-^Li2(-x)X - ^(x - ?/)(l - xy)X3 + ^Li3(-x) + 
+ ^(11 + 33x + 33x2 ^ 25x3)X - ^Cs - + 



+ — ^(18xy^ + 27x'y + 37x^ - 43y^)X^ + 



+ 



72y3 
77r2 



1 



37x 1307 



(20x^ + 25x + 11) + -^^(10 + X + 28x'^)X + — + 



216y2^ ' 36y2^ " ' ' ' 72y ' 864 

7 116 
+ 2%3 ~ + ^ ~ 3Li2(-x) - -XF + — TT^ + 

+ i(71x=^ + 138x2 + 147x + 43)X + ^(28^^ + x + 10) 



+ 



(A.17) 



E. 



[1] _ - 1) 



2 Li4(-y)-Li4 



Li4(-x) + (Li3(-x) - Cs)Yj - 

+ 



1 



y' L3 



^x3(Li3(-x) - Cs) - lLi2i-x)x^X 



+ 



24y2 



1 TT^ 

- ^(6a;^ + 3x2 - 7 - 3^^jj^3 _ ^4 _ _ 33.2 ^ 10x3)X 
-8(12x2 + 12yx - y2)(Li3(-y) - Ca + yLi2(-x)) + 



+ 



+ (3xy + 4y2 - 3x2)yx2 - 2(26xy - + 24x2)Xy2 



+ 



1 



, -10xy-27y^ + 36x2)^2 
Jbyx 

1 ^ 19x 3401 

— (5x + 9?/)F 4 h 

6yV ^ ^ 72y 2592 



36y- 



r(14xy-33y2 + llx2)X 



+ ivr 



-i3Li2(-x)(12xy2 _ y3 ^ ^3 ^ ^2x2y) + 
+ ^^^^^(Li3(-x) - Ca) + 3^(1 + 2xy)x2 + 



8y' 



+ 18^^'^'+1V^'""'^^'''^ 



6y 



ISy^x 



(6x2 ^ 22xy - 27y2)y 



^(x2 + 14y2 - lOxy + 12y)X + ^{^^^^ + 2x - 21) 



lx2-lx-^ + .Jlx-l 

24 54 27 V12 54 



[1] _ 



^X2 + 

722/ 72y 



(A.18) 
(A.19) 
(A.20) 



-33- 



For /i = 2 in eq. ( |2.23| ) and color factor tA'^'^ in eq. ( 2.30| ): 



= 4 



,{Ax - 1) (Lh{-x)Y - Lh{-x) - Li4 (^-^) ) 



1 

24 



1 

96 



(lOx^ + llx^ + 28x + 7)y^ + 



vr^ 1 

H (x^ - 2x + 2)xLi2(-x) + -(2x^ + 23x2 + 4x + 3)Li4(-y) 

6 4 

vr^ 1 
- 4^(22x^ + 75x2 _^ _^ -^5^y2 - -(x'^ + 9x2 _ 3^ _ i^^^^y - 



IT 



2880 



(9x^ - 182x2 + 113x + 24) 



^Y^U2{-x) + ^YLhi-y) + 



+ — 

yZ 



^(11x2 + 50x + 9) Uhi-y) + Li2(-x)y + ^xy2 



^(205x2 + 517x + 72)C3 



13(x - y) 



I8y 



Y^ + 



vr 



144y 



(445x + 208)y + 



288xy 



(1160x2 - 713x - 360)y2 + 



5x 



+ ^(3x + 2)^^ + 
1 



6y^ \ ' / ' 



2 , 2x + 135^^ ^ x(30377x + 30593) ^ 



3456y^ 



+ ITT 



1 



„ , (4x - l)xLi3(-x) - -(x'^ + 9x^ - 3x - 1)C3 - 
V 2 



x-y 
Ay 



+ -Xy2 + -^(11x2 + 50x + 9)Li2(-x) 
4 6y^ 

29x o 



13(x - y) 2 



6y 



yz ^ 



1 133 1 

+ 77^^' - T7^(-1160x2 + 713x + 360)y + — - — 
144y 144xy 54y 27 



(A.21) 



[2] 



2y3 



(4x^ + 12x2 - 3)Li4(-x) + (5x3 _^ 4^2 _ ^ _^ 7)Li4(-y) + 
+ (14x2 + 22x + l)Li4 ^-^^ - (5x2 ^ -^Q^ _ 2)yLi3(-x) 
- 2(2x2 + X - 4)xXLi3(-x) + ^(x2 - X + \)X^ + 



12 



+ y (2x3 - 6x2 _ y2,x + l)Li2(-x) + 

+ ^(x^ + 20x2 + 31x + 5)y^ - ^(lOx^ + 39x2 ^ ^ 6)Xy3 

7r2 7r2 

- Y^(5x3 + 10x2 + llx + 13)y2 - —(4x3 + 4x2 - 3x + 1)^2 - 

- (5x2 ^ _ ^^i^^^ ^ (^5^2 ^ _ 2)^gy _ 
^4 



TT 

1440 



(225x3 - 130x^ - 719x + 112) 



+ 



-34- 



2y2 



+ — -(liar - l)XLi3(-y) - {2x'^ + x + 2)X^U2{-x) - 

- 2(x - 2)XYU2(-x) - l{Ax^ + 7x + 9)X^Y + 

6 

1 TT^ 1 

+ -{x^-7x + 4)X2y2 _ (sx^ -5x + A)XY + 
+ ir. ( -i^^ + 4)rLi3(-y) - ^xY^U2i-x) + 7r2(y2 - X^) ) + 



+ 



1 

% 
1 

4y3 



(21x^ + 26x + 12)(XLi2(-x) - Li3(-x)) 



+ ^(7x2 + lOx + 4)X - ^(lOx^ + 5x + 1)X=^ + 
+ ^(20x2 + ISx^ - 12x - 4)X^Y - 

- ^ (443x^ + 706x2 ^ 299x - 216)C3 + 

(lOx + 3)(Li3(-y) + FLi2(-x)) - 
1 27r^ 

- -(lOx^ + X + i2)xy2 _ — (i9x + io)y 

6 9 



+ 



2y2 



15X + 28 o 



X 



18y 



144y' 



;(27x^ - 92x2 - 184x - 164)X'^ 



288y2 
1 



4(27x3 - 119x2 - 138x - 100)Xy - 7r2x(54x2 + 217x + 343) 



144x?/ 
71x - 24 



(27x-^ - 290x2 - 199x + 72)y2 - ^""^^ 



82/2 



X- 



2Ay 



-Y 



2592y: 



r(19139x + 18815) + 



+ in 



1 



L ( + '^^ + - 2)Li3(-a:) + Ca ) + 



+ 



1 

2?L 



(5x2 _2x + 5)Li3(-7/) + (5x2 + 3^ + 4)yLi2(-x) • 

- ^(2x - l)(x - 1)^3 + 1(7x2 + 7x + 12)Xy2 + 

7r2 7r2 1 

+ y(8x + 5)X - y(5x2 + 7x + 8)y + 

+ ^ ( 2xXLi2(-x) + ^(x + 3)y3 - l(y - 2)X''y\ + 



+ 



4y3 



(x2 + 6)Li2(-x) - ^(20x3 + 39x2 _^ _^ 20)^2 + 

+ ^(40x3 + 97x2 _^ _^ 38)Xy - 
7r2 

- — (x^ + 86x2 ^ -^53^ ^ gg^ 



-35- 



5(4a; + 7) . 



I2y 



96y 



-(73x2 + 74x + 100)X- 



^^(16x3 + 39a;2 + 3x - 8)y + -^imx'^ + 65a; - 24) 



(A.22) 



[2] _ 



1 



(2^2 + 4x + 1) (^u{-x) + Li4 (^-^^ + y Li2(-x)^ - 



- -(7^2 + 14x + 5)Li4(-y) + x(y - l)yLi3(-x) + 



TT 



24 



TT 



12 



+ ^(9x^ + l%x + 7)y^ - Cal^ -i-{x' + 2x- 1) f 



360^ ' 



+ 



+ ^ (^X{hh{-y) - Cs) + yLi3(-y) + \y^U2{-x) + ^X^F^ ) + 



+ 



x(3y + 1) 



4y3 



TT 



Li3(-x) - XLi2(-x) - -X'^Y + 

2 D 



-g^(3x2 + 4a; + 3)C3 + 9g^(18^' + 25x - 5)7r2 



+ ^^(2Li3(-y) + 2U2{-x)Y + XF^) 



(3a;2 + 12a;2 + 12a: + 4)X 



+ 



4yi 

X 

16y3 

-^(Gx^ + 12x2 - 15x - 8)F 
32xy ^ ' 

3x-4 



2 3x(y - 1) 

2 3x(y - 1) 
By 



By 



_ x(255x + 247) 

r2v + 



L2"7 (^^ " l)a^Li3(-x) - (3x3 ^ ^ 5^ ^ 2)C3 j + 
+ ^ (^4Li3(-2/) + Li2(-x)y - ^y + XF^ _ Jy3 ^ g^^^ ) ^ 



6 



1 



+ ^ (^-(% + 1)^^^ - + 2) + (7x" + 12x + 6)Li2(-x) ) - 



x(3y + l) 15X + B J_ . 
By3 ^+ IGxy ^By2^^ ^. 



(A.23) 



[2] _ X{X - 1) 



Li4(-x) - Li4(-y) + Li4( --j + — Li2(-x) + 



1 



1 



TT 



77r 



4n 



+ 7T7^ - + T^^ - ^(Li3(-ic) - Cs) - ^ 



24 



12 



360 



+ 



+ 



^(8x + 56)Li3(-2/) - 2(x + 7)Li2(-x)y + 



+ ^(17x + 65)C3-(a; + 7)xy2 



I 7(x-y) ^3 
72y 



-36- 




-37- 



For /i = 2 in eq. ( |2.23| ) and color factor Tr^^] in eq. ( 2.30| ): 



H. 



[3] 



1 



(5x3 _^ 5g^2 ^ _^ 13)Li4(-y) - -(lOx^ - 16 - 47x + 23x2)Li4(-x) + 



+ (16x3 + 7x2 _^ 323. _^ 8)Li4 



(12x3 + 2x2 + 28x + 7)yLi3(-x) + 



+ (12x3 + 21x2 ^ -^8^ ^ 5^ j^XLi3(-y) + -X'Y^ ) + 
+ (13x3 + 10x2 - 2 + 2x)XLi3(-x) - 



1 



(10x3 + 3x2 _ 9^ _ 8)x2Li2(-x) + 



+ 2(x - y)(2x2 + 3x + 2)XyLi2(-x) 
,2 



(44x3 + 109x2 + 68x + 22)Li2(-x) 



vr 
"6 

-(32x3 + 63x2 ^ gg^ ^ 32)Xy3 + -(x^ + 10x2 + llx + 4)x3y + 
6 6 
1 7r2 

— (19x3 + 16x2 + 41x + ll)y^ + —(10x3 + 11x2 + 13x + 8)^2 - 



— (17x3 + 30x2 _^ + Q)XY - —(9x3 _^ ^g^i _ -^3^ _ -^^^yj 
(14x + 13x2 + 11x3 + 5)C3^ + (11x3 _ 3.2 _^ 25x + 6)C3y + 



vr 

360 



(276x3 + 147x2 - 378x - 202) 



1 



x(2x-l) 4 
32y2 



+ 



^ ( (13x + 15)yLi3(-y) + ^(5x + 9)y2Li2(-x) + (2x + I)vr2(x2 - y2^ 
2y V 2 

1 



12y3 



(62x3 + 67x2 _^ g2x + 26)(Li3(-x) - XLi2(-x)) + 

+ (242x2 _^ 238x + 27)C3 + -(22x2 + 14x - 11)^3 

3 

- (21x3 + 48x2 + 39x + 16)x2y - 
2 



vr 



■(321x3 - 182x2 - 139x - 18)X 



+ 



+ 



12y2 



12 

(62x2 ^ 263x + 27)(Li3(-y) + yLi2(-x)) + 



+ (41x2 + 132x + 4)Xy2 _ ^(445x2 + 418x - 75)y 



+ 



18" 



72y 



: (172x2 +404X + 97)^2 + 



12y2 



(6x - \\)XY + 



+ 



-^(208x2 - 161x - 108)y2 _ -^^1 '^-2 
72xy ^ ' 72y2 

+ ^^(38x + 83)X + T-^(76x + 211)y + 



(18x2 - 108x - 49) + 



54y2 



108y- 



+ 



-38- 



( (x^ + 8x^ - 26x - 9)Li3(-x) - {x^ + 20x^ + 25x + 10)Li3(-y) + 

+ {3x^ - Sx^ -9x- 5)yLi2(-x) + 

+ ^{6x^ + 29x^ + 37x + 16)XLi2(-x) + 

+ ^{5x^ - _ -^5^ _ 8)jss:y2 _ 

- ^(lOx^ - Sx^ - 3x + 4)X^ - (14x2 - llx - 1)C3 - 

- — (3x^ - 15x2 - 33x -17)Y + — (Sx^ + lOx + 6)X j + 

+ ^ (^(5x + 3)y3 _ 15yX2y - 6(2x + 1) (xY^ + 7r2(X - F))^ + 



+ 



12y^ 



-(124x^ + 392x2 + 352x + 53)Li2(-x) + 



1 



+ -(64x^ - x^ - 38x - 6)X^ - 

- (62x3 + 117x2 + 60x + 13)Xy + 

7r2 

+ (124x3 + 29x2 - llOx - 77) 

64x + 25 1 

-F2 ___ (190^3 + 353x2 + 94x + 66)X 



24y 36y3 

r (190x3 + 116x2 - 335x - 108)r - 



36xy' 



121;r + 211 
108y2 



(A.28) 



r[3l _ 



1 



(x3 + 4x2 + 5x + 5)Li4(-y) - (3x2 + g^. _ i)^^^ _ 

- ^(2x3 + 3x2 + 93. _ 8)Li4(-x) - 3(y - l)xLi4 j - - 

2 \ y / 

- (2x2 + 4x + l)y(Li3(-x) - Cs) -ix + 3)j/2yLi3(-y) + 
+ (3x2 ^ _ i)xLi3(-y) + ^y^Y^U2{-x) + 

+ (x3 + 4x2 + 8x - 2)XLi3(-x) - 2x(y - l)Li2(-x)Xr + 

7r2 1 
+ —(4x3 - x2 - 14x + 2)Li2(-x) - — (x3 - 3x + 1)F^ - 



24^ 

^(2x2 + 7x + ll)X2Li2(-x) + ^(5x2 + lOx - 3)X2y2 



1 7r2 
- -(x3 + 3x2 ^ 3^ _ i)x3y - — (y - l)(x2 - 3x + 1)XY - 

D o 



(2x3 + a;2 + 2x + 15)^2 + — (x3 - ^2 - 5x - 4)y2 _ 

12 4 

- ;^(2x3 + 8x2 ^ ^ - ^(2x - 1)X^ - 

6 48 



TT 

360 



(20x3 - 85x2 _ 248a; + 20) 



+ 



-39- 



+ 



4y3 



+ 



-{x^ + 4x- 6)(Li3(-x) - XU2i-x)) + -x'X-" + 

1 TT^ 

+ -(9x + 5)X2y + 9(x - y)C3 - -^(14x + 47)X + 
i(7x + 9)(Li3(-y) + U2{-x)Y) + l(10x + U)XY^ + 



+ 



TT 



+ —(2x3 + 12^2 + 9x - 2)y2 - — ^ 
8xy 2%^ 



(Gx^ + ISx^ + 14x - 9) 



-^(23x + 26)X - -^(46x + 47)F + 



+ ivr 



TT 



^ (^-(x^ + 2x2 + X + 4)Li3(-y) - '—{x + 3)(x - y)X + 

+ (x=^ + 2x2 + 4x - 3)Li3(-x) - (x^ + 5x - 2)C3 + 

+ ^(2x3 + x^ - x + 6)X3 + i(y - l)(x - ifXY'^ 

- I (2x2 + 3x + 3)XLi2(-x) - 



(x^ + x2 - X + l)yLi2(-x) + 



TT 



+ y(x''-2x^-7x-2)y 



X - 2, 
T2^ 



+ ^ (^-3(2x^ + 4x + 5)Li2(-x) + -(3x^ + 6x + 2%)X^ + 

1 TT^ \ 

+ -(x2 + 20x + - _(2x - 1) j + 

+ 7r^y2 _ -l-(24x3 + 25x2 + 4x + 6)X - 

- —^(24x3 + 38x2 ^ -^5^ _ ^ -^(41x + 47) 
12x?/^ 12y^ 



(A.29) 



J, 



[3] _ {X - 1)X 



hu{-x) + Li4 (^-^^ - Li4(-y) - yLi3(-x) + 



+ yLl2(-x) 



TT 



4n 



xy^ + — + — y^ + 



24 



12 



77r' 
360 



+ 



+ 



+ 



24y^ 



(GX^ + 15x2 + 15X + 5)X2y + + 8x2 ^ -^9^ ^ -^g^jsj. 

-(x2 + 6x - l)(Li3(-y) + yLi2(-x)) + ;^(2x - l)xX^ - 



3yH 



+ 



- -(x + 9)(2x - l)Xy2 + (5x - 1)C3 



+ 



+ ^{{y- l)(Li3(-^) - ^Li2(-x)) - ^(5x - 21)y ) - 



-40- 



- -Y^ + ^^(6x3 + 41x2 + 48x + 19)X2 + 
9 24y3 

+ ^77^(363;^ + 108x2 + 37x - 29)Xy + ~ 



72y2 ^ ^ I8y2 

1 _ 9 ^„ ,_._2 (20X-19) 



+ — (6x'' - 5x^ + 13x + 12)y^ + 

24xy ^ ^ 

2 

36y2 



72y 



TT 



(9x2 ^ 27x2 + 14x + 11) + 
^^^(Li3(-x)-C3)-- 



1 
3y 



2(2x2 + 9x + l)Li2(-x) + 



+ (2 (4x2 + 6x + 3)xXy - (6x2 + 5^2 ^ ^ ^ 3^jsj.2~J ^ 



24y 

+ (102x2 + 143x2 + 106x + 29)X + 

+ -^(102x2 - 11x2 - 179x - 72)y + 

72xy^ 

X — n tt'^ , . 17x — 19 



(A.30) 



2 ^Li4(-x) + Li4 ^-^^ - Li4(-y)^ (x - l)x + 



vr 



+ 2(-Li3(-x) + C3){x - l)xy + yla; - l)xLi2(-x) - 

- J(x - i)xy2 + :^(x - i)y^ + ^(x - i)xy2 - ^(x - i)x 



12 



180 



+ 



+ 



1 

%2 



1 



-x2Li3(-x) + x2XLi2(-x) - -X^Y - 



- (15x2 + 15x + 1)C3 + ^(10a;2 + 3)X 

-^(x2 + 14x + l)(Li3(-y) + yLi2(-x)) - 

1 7r2 

(8x2 ^ ^ 7)xy2 + (4y + i)y - 



24y2 

1 
8? 



12y 



vr 



(4x^ + X - l)Xy - —(18x2 + 43x - 9) - 



,2 8x + 29 



y--^(2x + 7)X + 



(x — l)x 



(-Li3(-x) + Ca) + 



+ ( (2x2 + 15x2 + 15x + l)Li2(-x) - ^(4x2 + 1)X2 + 
oy \ 8 

+ -(4x2 + 11x2 + lOx + 2)Xy - — x2^ - ly±ly2 ^ 
4^ ^ 6 y 24y 



- 41 - 



+ ^(5-^ + 8x - 1)X - ^(21.^ + 35x + 12)Y + ^ 

72y 36y 36 54y 24y 54 

+ ^(^(5 + 4x)y + (l-4x)X-- 

3 in eq. ( |2.24| ) and color factor Tr'^^ in eq. ( 2.31j ): 

^ ' -3x^ - 31x2 _2x- 4)X2 + 



48xy2 

1 , 9 9x 175y 25 

+ 55xy - 16y^ + 60x^)X + ^ + 

A8xy^ ^ ^ ' 108x 864 

+ mi —^{-3x^ - 31x2 - 2x - 4)X + ^— (55xy - 16y2 + 60x2 
\24xy^ 48xy 

^(4x3 + + 7x2y + 3xy2)x2 - 



16y 



^ r(6x2 + 6y2 + llxy)y2 _ 1(1 _ y)YX 



16x2' ^ 8 

2 , , 3(1 -w) 7r2^ X 24 139 



7r/5x + l 7w + l 1,9 9 

+ i- ^X - -^^-^Y + — 18y2 + 31xy + 12x2 

y x^ Axy 



y .(6^2 + 62/2 + iixy)y2 + y-^YX - %r^X + 



16x2' 8 16y 

+ ^^(6x2y + 4xy2 -y^ + 2x^)X'^ + (y - 1) f —Y -—)- — + 
16y2' ^ ^ ^ ^ '\Sx 16/32 



— (1 + 3y) + —^(3x2 + 2xy + 2y2)X2 + 



I . J_(y2 + + 2x2)X + ^-±^Y - -^(-12y2 + 4x2 
y^ x^ 2xy 

5x2 + x + 2^^2 x2 + 23x + 10,, 397y 143 
24xy2 24xy 216x 108 

/5x2 + x + 2,^ x2 + 23x + 10 

+ ITT ^ X H 

V 12xy2 24yx 
^18a_4^^__^^^^,^g^,^^^^^^,^ 
4' 4x 72x 

12x 54x 12x ' 
0. 



3 in eq. ( |2.24| ) and color factor tA'^'^ in eq. ( 2.31| ): 



1 

xy 



- — (SOrc^ + 39x^ + lOx + 4) (^{X - Yf + vr^ 

- -^(eOx^ + 55x - 16)X + ^(180x2 + Ulx + AO)Y + 
48 144 



OK 

+ — (x - 56) 
864^ ' 



ITT- 



- 2) 
36xy 



6 — a; 
16x2 



y^ + i^x^ + ^(4x-i)((x-y)^ + .^) 



6x + 11 1 2 115x - 24 



TT / X — 6 



1 

16 



x^ 
X — 6 



y-1 



x + 



y2 



y-1 



13x - 18 
3xy 

1 



x^ 



^4x + l^ 

y 

.IT f X — 6 



1 

xy 



X^ + -(2x^ + 4x + 1) (X - y)^ + vr^ + 



(4x2 _ 5a; _ i2)y + 



17 

2y. 



-^8 
1 

24xy 
1 

~ 24e 



y + ^x + 3^ 



xy 



(12x2 + 3x - 10)X - ^(36x2 + 49x + 22)y + ^(286x + 397) 
3 9 



(6x2 + 3x + 2)((X-y)2 + ^' 



l Ox + 13 

18xy 



18x + 35 



y-1 



x^ 



y2 + (3x + 2)((X-y)2 + 7r 



3y 



X + — (12 + 19x + 6x2)y + ^±}^ 
xy 9xy 



TT X 



— I ■ 



-X + 



12x 



y-1 



10 



y 



llx + 18 



3xy 



TT 



3 in eq. ( |2.24| ) and color factor Tr^'^^ in eq. ( 2.31| ): 



y-1 
8x^2/ 



(1 - xy)y2 



48y2x 



(6x3 ^ 25x2 - 7x + 4)X2 - ^^"^ + 



+ ITT 



^^(30x3 + 51x2 + 19x + 4) ( (X - y)2 + 7r2 ) + 



18x 

7(y - 1) 



y 



6y 



\ 4x^7/ 



(1 - xy)y 



^ (6x3 + 25x2 - 7x + 4)X + + 



^(x2-2x-2)y2 
8x2y ^ ^ 



24y2x 
X — 1 



18xy 



((X - y) 



2 + 7r2 ) + 



J. 



[3] 



Kl 



[3] 



L 



[3] 



A 



[1] 



16y 



12y 



12x 



X — V 



24:xy 



/ 9 39a; + 94 



4x2 



72x 



12xy 



+ ivr 



X - 2 



24y 

2x2 
-X^ + 



y + 
y-i 



12xy 
y2_ 



r(5x^ + x + 2)X 



55x + 94 
72xy 



1 



8y2 8x2 
19x + 36^^ 19(y-l) 



+ 

+ ITT 
1 

36x 



24x 

y-1 



-Y 



2iy 



(3x + 2)((X-y)2 + ^2'j ^ 



4x2 4y2 24xy / 

(y + ivr) . 



For /i = 4 in eq. ( |2.25| ) and color factor Tr'^^ in eq. ( 2.3l| ) 



1 



^(x^ + 3x2 + 3a; - 2) (jt'^X^ - Li4(-x)) - XLi3(-x) + 

+ (1 - xy)X^ Q(2 + x)Li2(-x) + ^{y- l)X^^ + 

1 13 

+ -(15X + 5x3 ^ ^5^2 ^ j^^^j^ ^ —TT^X^ + 

2 D 



TT 



+ -^i^x + 1)(19 + 24x + 9x2) 

7r2 7r2 1 

+ — Li2(-x) + —YX YX^ + 

6 ^ ^ 6 12 



+ 



+ 



1 



288y3 



4(50 + 108x + 126x2 + A5x^)X''^ - 
- (579x + 795x2 _^ 297x3 + 265)7r2x 



+ 



+ 



48y2 



(44 + 73x + 35x2)(2XLi2(-x) + yx2 - 2Li3(-x)) + 



+ (241x2 + 470x + 205)C3 



+ 



288r 
1 

288y2 
7r2 



■(1530x2 + 505x3 _^ _^ 9i8x)X^ 
(18x2 - 431x + 501)X + 



+ ^:3-^(18 - 73x- 45x2) 
iZy 



475 36077 



288y 3456 



+ 



(A.48) 



(A.49) 



(A.50) 
(A.51) 



-44- 



TT 



4 [y 



^ (2 + x)(l - xy)XU2i-x) - (1 - xy){l - y)X^ 



- 2Li3(-a;) + —X + {bx^ + ISx^ + 15a; + 7)C3 



+ ^(50 + 108x + 126x^ + Abx^)X'^ + 



1 



+ g^(44 + 73x + 35x2)(Li2(-x) + XF) - 
-(505x^ + 1530x2 + 368 + 918x)X - 



36y< 



24y- 



r(21x2 + 50a;+45) + 



475 467 1 



+ 



+ 



36y2 72y 4j ' 



(A.52) 



1 

4y3 



[y - l)(4a;2 - 5x - 5)Li4(-x) + ^(a;^ - 3x - 3)X^ - 

- 2(4a; - 4^^ - + l)XU^{-x) - 

- — Li2(-x)(25x^ + 44x2 + 31x + 14) - 

O 

- (llx^ + 17x2 + 13x + 9)C3X - y(4x^ + 3x2 _ g^. _ i2)x2 + 



+ ^(653x^ - 133x2 - 1241x - 343) 



+ 



+ 



1 

2? 



X 



-x(x + 13)Li4(-y) - (7x^ - x + A)hu + 
+ (3x + 4 + 5x2)(Li3(-x) - Cz)Y - 2(2 - X + 3x2)XLi3(-y) 
+ ^(6x2 - 5x - 2)X2Li2(-x) - 2Li2(-x)(2x2 + x + 2)XY - 



- 2yY'^U2{-x) + 4y2yLi3(-y) - -(-3x + 13x2 + 8)X2y2 + 

+ J (5x2 ^ 4^ ^ 5^y3jsj^ - i(l + ^x)YX^ - 

o u 

1 7r2 

- ^(5a;2 - 5x + 2)F^ - y (8x2 + + 7)^2 ^ 



TT 



+ :-(6 + 4x + I3x2)yx 



+ 



+ — 

y6 



^(4 + lOx + 3x2)(XLi2(-x) - Li3(-x)) - 
- ^(6x + 32x2 + lh)X^ + 
+ Y^(1509x + 443x^ + 443 + 1617x2)^3 + 



+ 



TT 

288 



(^477 + 891x + 999x2 ^ 7333.3 



+ 



+ (x(llx - 19)yx2 - 2(x2 + 9x + 10)7r2y) - 



-45- 



- [3xiUs{-y) + YU2i-x)) + — (27x + 20)XY' 
^Y^ - 3^^(246a;2 - 215x^ + 246a; + 1U)X^ - 

^ (4x(145 + 91x)XY - (311x2 - 340x - 455)7r2) + 



288y- 

+ 1^^^^^"^ + " 21^ ^^^^"^^ + ^^^^"^ + ^^^^^ + 



3(2: -y ,. 139 14135 

+ —Y h + 

8y 72y 2592 



+ ivr 



^(x^ - 4x2 - llx - 5)Li3(-x) - ^(-7x2 + 2x^ - 6 - 14x)X + 

+ \{y^ + 2x2 + 4x)C3 - ^xixy - 2)X^^ + 
+ ^ (-(^ - 5)^Li3(-y) + J(2x + 5)(y - l)XLi2(-x) - 

- x(2x - l)yLi2(-x) - Y^>"(17x + 3x2 ^ _ 
-^(l-2x + 3x2)Xy2^ + 

+ - (-^xY^ - IxXY^ + + 1)X^Y - ^xX 
y\ 12 2 4^ ' 2 

2 

- -^(9x2 + 14x + 8)Li2(-x) + -^(23x^ + 93x2 _^ ^g^^ _^ _ 
Ay^ Uoy^ 

- ^(13x2 + 7x + 2)X2 + ^—(29x2 + 39x + 40)XY + 
6y^ 12y^ 

21x + 10^2 , 1 ,.n«™3 ...2 

+ ^^^^(100^^ + 264x2 + 27 + 245x)y - 
^ (431x2 + 2209x + 863) 



+ — y2 ^ ^^(306a;3 - 10x2 - lOlx - 144)X + 

I2y 72y^ 



216y 



x(x2 - 4x + 1) \Lu{-y) - Li4 (^-- j + (Li3(-x) - Ca)^^ - 

24 12 6 7 

+ ^x(5x2 + 8x - 2)Li4(-x) - 2x'^XLi3(-x) + 

3 2 3 

+ ^Li2(-x)X2 ^ !L^(^2 ^ 4^ _ l)Li2(-x) - - 

7r2 



(A.53) 



(x'^ + 9x2 ^ 9^ ^ 3)^^^ ^ l2^^ 

(16x^ - 5x2 - 40x - 11) ^ 

360^ 



-46- 



+ 



9x{5x'^ + 1 + Ax){U2{-x)X - Li3(-x)) + 
+ 3(15 + 29x + 21x^ + 29x2)C3 + x{l + 4x + 9x'^)X^ + 



+ (3 + lOx + 13x^ + 9x^)7r^X 



+ 



+ 



4y2 



3x{x - l)(2Li3(-y) + 2U2{-x)Y + XY^) - 
1 



12 



x{3x + l)y(3X^ + 47r^) 



^(17x^ + 6 + 34.)X^- — 



{2x + llx^ + 3)y2 - 

- ^^(l2x{5 + lla;)Xy - (51a;^ + lOOx + 29)7r2) + 

1.2 3(x-u),, 17 187 



^ (^6x(l - 4x - a;2)(Li3(-x) - Cs) + 



+ 



+ 



1 



24y3 

8^ 



+ 3x'^U2i-x)X + X'^TT^X - x'^X'^ + oCs + 



9x(a;^ + 4a; + 5)Li2(-a;) + Qx'^X^ - (3a;^ - - 8x - 3) 
16y3 ' 



(3x + l)Xy + -^(5x2 + 4 - 2x)X + 
Iby^ 

1 o _ 1 



(8^2 + 5aj + 3)F + ^(3x + l)(14x - 9) 



1 



2(9x^ + 27a; + 27x^ + 14)X^ - (47 + 81a; + 27a;^ + 81a;^)7r^X 
17 . 1 



+ §(Li3(-a;) - U2i-x)X) - ^C3 - ^YX' + 

1 -jr2 

+ ^;r^(43a;^ + 90^ + ISSx^ + 35)X^ + ^(77^^ + 139a; + 32) + 

72y3^ 216y2^ ' 



+ 



1 



216y2 
+ in 



(53 + 400x + 125x^)X + — + 



37 1307 



72y 864 



+ 



\u2{-x) - \XY - ^{9x'^ + 27x + 27x2 + M)^^ + 



+ ^TT^ + ^(43x3 + 90x + 135x2 + 35)X + 
16 36^* 



+ 



1 



216y- 



r (125x2 +400X + 53) 



x(x — 1) 

y3 



2 Li4(-x) + Li4(-y) - Li4 



+ 



TT 



+ 2{X F)(Li3(-x) - Cs) - yLi2(-x) 



Ixy3 + ly4 1^4 !Ly2 + l;r' 

3 12 12 6 3 60 



+ 



[1] 



[1] 



+ 



144yJ 



48(15x + 15x2 + + 2)(Li3(-x) - Li2(-x)X) + 
+ 2x(13x2 - 24 - 24x)X^ - 2{x'^ + 3x2 + 25 + 3x)C3 



37r2(13x^ + 9 + Tlx + Tlx^)^ 



+ 



+ 



6y2 



-2(x2 + 14x + l)(Li3(-y) + yLi2(-x)) + -(Gx^ + 17x + S)YX^ 



(2x2 + 17X + 3)Xy2 _ y (2x2 _ _ j^y 



+ 



+ _Ly3 
12 



^(123 + 312x + 244x^ + 54x3)X2 + 



72y 



+ 



TT 1 

(lOSx^ + 287x2 + 128x + 33) (44x + 27x2 - 19)y2 _ 

36y 



144y2 

X 



1 



+ 



18y2 
9x + 4 
6y 



(43 + 76x + 27x2)Xy + 7tt^(353x^ + 478x + 191)X + 



216y2 



+ in 



19 3401 

Y H h 

72y 2592 

„3 , c;i™2 



f ^(21x^ + 51x2 + 39x + 17) - Li2(-x) + 
6y^ \48 



X 



+ -(13x + llx^ + ll)X^ ) + 
^ x(2x + 5) _ (x-l) ^^2 



12y2 



12y 



36y 



r(38x2 + 37 + 74x)X 



1 (18x + 5x2 + 19)y - + — 

^ ' 36y2 216 



18y2 
2 5 



vr 



-r^ - — -X + ivr —X 



24 



27 54 



12 



54 



^{xYX + y2y + 7r2) + i ^{xX + {y- l)Y) . 

For /i = 4 in eq. ( |2.25| ) and color factor Trf2] in eq. ( 2.31| ): 
1 



A 



[2] 



4y3 



-4(x - 4)x(Li4(-x) + Li4(-y)) 



(3x^ + 4x2 ^ 23x + 2) Li4 



"6 



Li2(-x) 



- 2(x^ + x2 + llx + 1){X - y)Li3(-x) + 

+ — (3x^ + 14x2 + 7x + 6)X^ + -(x^ + 5x2 + X + 2)y4 
12 6 



(x^ + lOx^ - 7x + 4) XY^^ - TT^X^ + 



IT 



+ ^(x - l)(x2 + 9x - 2)y2 - 4(6x + iXsY 



(A.56) 
(A.57) 
(A.58) 



-48- 



- (9a;^ -6x + 5)C3^ - ^(12^;^ " ^^^^ + l^^x - 21) 



4y V2 



^{X -Y)^U2{-x) + ^X^Y^ 



- i (^l(3x + 5)X^Y + ^XY^ + i(X - F)Li3(-y) - 



1 



- J (16x2 + Tlx + 25)X^ + J (24x2 + 95x + 41)Xy2 _ 



vr 
24 



(176x2 + 683x + 267)y - 



^ 7r2(40x2 + 203x + 67)X - ^yCs 



117 



+ 



48 

3y V 6 ; ^ 288y2 



+ 



(SGOx^ + 1073x2 + lOOlx + 528) - 



( x(360x + 713) (X - y)2 - 368^2 + 560Xy - 1160^2 ) + 

288y \ / 

1 / 1 , , 1 , , \ 30593x + 31 



30593X + 30377 



18y ^ ' 86Ay 



, Il7r2 

Cs H h 

4^^^^ 288y 



(A.59) 



[2] _ 



1 



(3x3 _ -^22; _ 4)Li4(-x) - (x2 + 22x + 14)xLi4(-y) - 

- (7x3 _ 3.2 ^ 4^ ^ 5)Li^^_^^ _ (7a.3 + 3.2 _ 4^. ^ l)XLi3(-x) + 



+ (6x3 _^ 3^2 + 9^ + 5)yLi3(-x) + y (23x3 + 39x2 _^ 24x + 13)Li2(-x) + 

+ ^(x2 + 4x + 2)X^ - ^(3x + 2)(x2 - 4x + 2)y^ + 

1 7r2 
+ -(8x3 ^ ;lq^2 ^ 23x + 14)Xy3 -—(x^- Yix^ - 4x + 12)^2 - 



TT 

12 



(15x3 + 27x2 + 36x + 17)y2 - (6x3 _^ 3^2 + 93,.^ 5^^^^ _^ 



+ -(12x3 + 7x2 + 83. + 11)^3^ 



TT 



1440 



(592x3 - 551x2 - 1570x - 343) 



+ 



+ 



2y2 



1 



(5x^ - 2x + 5)XLi3(-y) - -(4x^ - x + l)X^Li2(-x) + 
+ (4x2 + 3x + 5)XyLi2(-x) + ^(5x - l)x3y + 



-49- 



1 TT^ 

+ ^(lOx^ -3x+ 11)X^Y'^ - y(7x2 + X + 4:)XY 
+ ^ ({4x + 5)YUs{-y) + h^^U2i-x)j + 
x{6x^ + l)(Li3(-x) - XU2{-x)) - 



+ 



1 

4y3 



1 



:^(14a;2 + 40x + 29)^3 + J (4^2 - 3)^2^ 
18 2 

75-2 

— (176x3 + 853x2 _^ 7g2^ _^ 477^^^ _ 
^(1235x2 + 1426x + 443)C3 + 

oD 



-(3x + 10)x(Li3(-y) + Li2(-x)y) - 



1 TT^ 

- -(33x2 ^ ^ 20)Xr2 + y (lOx^ + 28x + 15)y 
^ (72x'' + 145x3 - 304x2 - 422x - 144)X2 - 



14Ay 

_ 28x + 15^3 _ 1 _ _ _ ^Q^xY + 

18y 8y2 

7^2 1 
+ —^(144x3 + 146x2 - 83x + 167) + —^(216x2 + 917x + 782)X - 

1 o 9 24x-71_ 18815X + 19139 



-(72x3 - 199x^ - 290x + 27)Y^ + — — Y + 



144xy' ' 24y 2592y2 

+ m 



+ TT^ ( (x - ll)xLi3(-y) + 2(2x - l)xyLi2(-x) + xX3 + 

7r2 7r2 \ 

+ (3x2 _ + i)xy2 + (3a; + 2)X - -7-(3x2 - llx - 8)Y ) + 
+ 1 (-2XLi2(-x) - Jx2y + ixXF2 + ^y3 + _ + j + 

+ ^ (^(12x2 + 26x + 21)xLi2(-x) - |(20x2 + 40x + 41)^2 + 

+ ^(60x3 + 136x2 + 119x + 40)Xy - 



7r2 \ 
— (112x3 + 79 + 258x + 263x2) j 



^%i^y2 _ ^(100x2 + 74x + 73)X 
6y 72y3 ^ ^ 

1 , o 1 



(100x3 + 138x2 + 119x - 27)y + — — ^(1340x + 1421) 



72xy2 ^ ^ 216y 



(A.60) 



-50- 



c, 



[2] _ 



1 

2? 



-(Sx^ + Mx + T) Li4 



^Li,(-.)-lxy3)_ 



- (x^ + 4x + 2)x{hu{-y) + Li4(-x) + yLi3(-x)) + 
+ (2x^ + 6a; + 3)xXLi3(-x) - ^(2x^ + 2a; + + 



+ — (a;2 + 4x + 2)F^ + 24 ^(^^c^ + lOx + 5)^^ + 



+ 3^a;(a;2 + 6x + S)^^ + (x^ + 5x + 3)C3^ + 



+ (x - y)xC3i^ + T^(41x^ + 113x2 +40x - 11) 
+ ^ (2^ - Y)U^{-y) -\{X- YfU2{-x) + 



+ 



1 3 vr^ ■ 

+ -x^F + -x^y^ xy 



+ 



y.i 



--(6x2 ^ ;l2x + 7)(Li3(-x) - XLi2(-x)) + 



^ „2 I q™ I 1 \ -1^3 /I i™2 



+ — (6x^ + 8x + 1)X^ - -(llx^ + 24x + 15)C3 + 
z4 8 
2 

+ ^(2x + 3)X2y + J-(12x^ + 19x2 + 5x - 3)X 
8 24 



+ 



+ 



+ 



x(3.T + 2) 



2y2 



TT 



Li3(-y) + yLi2(-x) -'■jY + -xy2 



32y3 
2 



?x^ + 31x2 + 42x + 21)X 



2 x(8x + 15) 



XY 



96y2 

1 



(24x'^ + 69x^ + 86x + 29) + 



(16x2 + 31x + 3)X 



32xy 

4x-3„ 247X + 255 



16y 

;8x^ + 15x^ - 12x - 6)y' 



8y 



-y 



128y2 



+ 



|.(u,(-.) + Li3(-,)-ix» + ix=r-^x) + i-<3 + 



+ 
+ 



x(2x + 3) 



(6Li2(-x) + 6Xy - 3^2 - 7r2) 



24y3 

x(2x + 3) ^ 3(x - y) 



+ 



8y3 



8xy 



y 



32y- 



r(9x + 5) 



(A.61) 



D 



[2] _ X(X - 1) 



Li4(-x) + Li4(-y)-Li4(^-- ) +X(-Li3(-x) + C3) + 

.2 



1^" ^ . , 1 ,^4 1 3 



+ y(Li3(-x) - C3) + y Li2(-x) - -X^ + -XY' - 



TT 



TT 



^Y" - -rX^ - ^Y' 



7 



:7r 



+ 



+ 



3y2 



24 6 12" 120 

(7x + l)(Li3(-x) + Usi-y) + (y - X)Li2(-x)) - 



- 51 - 



- -^(5ic^ + 43x + U)X^ + ^(3a;2 + 37x + 10)XY^ 



2 

- ||(7x2 + 76x + 21)y - ^(lOx^ + 143x + 37)X 



TT 

48' 



+ 



TT 



216y2 
1 



(lOSx^ + 321x^ + 332x + 191) + 



+ — |^x(36z + 71) (X - Yf - 35^2 + 18XY - 71Y^ ] + 

216x + 53 9x - 7,^ 863 
+ X Y + — — + 



+ ITT 



E, 



[2] 



6y 
x{x — 1) 



216y 9y " ■ 864y 

^-y((v ]r^2, 2\ 13 31 

V^-^^ +^ J -144^-1%^- l%^+216y. 



221 



-2 (^Li4(-x) - Li4 (^-^ j + Li4(-y) j - y Li2(-x) + 
+ 2{X - Y){U,{-x) - Cs) + ^{X' + Y') - 



3 6 ^ 60 



+ 



+ 



3y3 



TT 



x^ + 15x2 + 15x + 2)(Li3(-x) - XLi2(-x)) + 

2 



+ — (42x^ - lOQx^ - 116x + 27)X + 
48 

+ ^(13^2 - 14x - 19)X3 + ^(^2 + 2x - 23)C3 



+ 



+ 



3y2 



-(x^ + 14x + l)(Li3(-y) + yLi2(-x)) + 

+ -^(llx^ + 16x + 8)x2y - i(5x2 + 57x + A)XY'^ + 



(A.62) 



+ — (13x2 + lllx + 2)y 



+ 



13X + 6 3 
72y 



72y 



:(54x^ + 235x2 + 338x + 139)^^ 



+ 



+ 



36y2 

144y2 
1 

21%2 
1 



(54x2 + 191x + 125)Xy 



72y 



(54x2 + 147x + 2)y2 + 



(108x^ + 362x2 + 117x - 33) + 



3x + 2 



2y 



y + 



(324x2 - 47x - 191)X - 



4085 
259% 



+ 



3y' 



3 , -Li2(-X) + -(6x2 ^ ^ ;^7^j^2 ^ 



+ !l!(48x3 + 135x2 + 126x + 47)^ + I^y2 + 
48 J I2y 



-52- 



+ -^{Wx'^ + 15x + 8)XY + —^(5x2 - llx - 7)X + 



1 



+ 7^(5x2 + 12x + l)y 



587x + 407 
216y2 



[2] 



vr 



[2] 



[3] 



-G 



9y 

[1] 



+ i- (x + y)-- , 



(A.63) 

(A.64) 
(A.65) 



For /i = 4 in eq. ( |2.25| ) and color factor Tr'^^ in eq. ( 2.3l| ): 



1 

2^ 



1 



(16x^ + 47x2 - 23x - 10)Li4(-x) 



- (13x3 + 31x2 + 56x + 5)Li4 ^-^^ - (Sx^ + 32x2 + 7x + 16)Li4(-y) - 

- (lOx^ + 25x2 + 20x + l)XLi3(-y) - (x^ - x2 + 28x + 2)XLi3(-x) + 
+ (Bx^ + 15x2 + 39x + l)yLi3(-x) - {5x'' + 9x2 _^ 3^ _ 3)xyLi2(-x) - 

- ^(2x + 3)(3x2 + 5x + 4)x2Li2(-x) - ^(4a;2 - llx + 11)X^ + 

+ —(25x3 + 53x2 _^ _^ ll)Li2(-x) - 
6 

- ^(13x3 + 21x2 + 3x - 13)x2y2 _ 1(9^3 _^ 33^2 _^ _^ 22)X''^Y + 
+ ^(34x3 + 86x2 _^ ^Q3^ _^ 18)XY^ - ^(lOx^ + 22x2 + 47x + 2)y^ + 



+ — (x + 3) (14x2 - 13x + 4)^2 + —(4x3 + 8x2 _^ 4^ _ 5)^^ 



12 



(7x3 _^ 9^2 _^ 3Q^ _ 5^y2 _ ^g^3 _^ ^g^2 _^ 42x + 2)C3y + 



+ (9x3 _^ ^4^2 _^ 34^ _^ ^^^^^ _^ —(10x3 + 334x2 _ j2ix + 184) 

360 



+ 



1 



+ ^ (^2(15x + 13)yLi3(-y) + (9x + 5)y'Li2(-x) + 2(1 - y){X'' - Y'')7t'' ) + 
(53x3 _^ 352^2 _^ 3g2^ _^ 124)(Li3(-x) - XLi2(-x)) - 

12y3 _ 

- (19x3 + 43x2 + 47x + 19)x2y + 
+ ^(68x2 - 7x - 35)X3 - (52x2 + 83x + 62)C3 - 



+ 



12 



(585x3 _^ 2038x2 + 2147x + 780)X 



+ 



+ 



12y2 



-(27x2 + 263x + 62)(Li3(-y) + yLi2(-x)) 



(4x2 + 132x + 41)Xy2 - ^(309x2 - 646x - 307)y 



+ 



-53- 



+ —Y^ ^(108x^ + 949a;2 + 1193x + 487)X^ 

18 72y3 

^ :{3Qx^ + 185x2 + 108x + 10)XY + 



122/2 

+ ^(108x^ + 450a;2 + 363x + 250) + (211x + 76)Y 

727/^ 108y 



1 x 
— (108x2 + 161x - 208)y2 + 
72y ^ ^ 108y2 



(108x^ + 161x - 208)y^ + :^^^(211x + 121)X + 

^ (^{7x^ + 16x2 + llx - l)Li3(-x) + 

+ (Sx''^ + 18x2 + 21x + 12)Li3(-y) - 

- ^(16x^ + 37x2 _^ 29x + 6)XLi2(-x) 

- 2{y - l)(2x2 + 3x + 2)yLi2(-x) - 
7r2 

- — (8x^ + 28x2 + 32x + 13)r + 
+ ^(9x^ + 36x2 _^ _^ 22)Xy2 _ 

- ^(6x^ + 9x2 + 9^ + g^x^ - (4x2 + 8x + 1)C3 ) - 
- i^(2x2 + 3x + 2)X - i Q(7x + 11)x2f + 

^ ( (26x3 + 62x2 + 67x + 62)Li2(-x) - 



12y^ 



- ^(151x3 + 252x2 + 228x + 86)^2 + 
+ (57x3 _^ jQ4^2 _^ ^ 18)Xy - 

7r2,_ o „ n A 21X + 8, n 



36y 



- —(44x3 + 17x2 - 80x + 9) Y^ 

12 / 8y 

^ (286x3 + 314x2 + 133x - 30)X - 



,3 ' 



1,9 , 83X + 38 

r(286x2 + 371x + 238)y- 



36y2^ ' 54y2 



.[3] ^ J_ 
4 2y3 



(A.66) 



-3(x - y)xU4{-y) - ^{8x^ - 9x2 _ 3^ _ 2)Li4(-x) - 

- (5x3 ^ 5^2 + 4^ + i)Li^ _ (4^3 + ^2 ^ 2x + l)XLi3(-y) - 

- (x3 + 5x2 + 4^ + 2)XLi3(-x) + (2x3 + 6x2 + + ij^-gX + 

+ (x - y){l - xy)F(Li3(-x) - Ca) + ^(2x3 + ^2 + 5x + 4)X2Li2(-x) - 

- (x3 - x2 + X + l)XyLi2(-x) + 

7r2 1 
+ y (x3 - ar2 - 2x - l)Li2(-x) + ^i^x - l)(x - y)xX^ - 



-54- 



24 3 

- ^(4x^ - + 2x + 3)X^Y^ + l{7x^ + 12x'^ + 12a; + 4.)XY^ 
4 6 

- j7:{9x'^ + 8x + l)xX^ + —(4x2 - 2x - l)xXY - 



- — (4x-^ + + 7x + 3)y^ + ^(62x^ + 140x + 79)x 



+ 



4y3 



X 
~~2 



4y 



3(5x2 + 4x + 2)(Li3(-x) - XLi2(-x)) - 

- ^(8x2 + 7x + 2)X^ + J(6x2 + llx + l)x2y - 

- y (22x2 + 28x + ll)X + (2x - 7)C3 



(9x + 7)(Li3(-2/) + yLi2(-x)) + 
+ J(13x + 10)Xy2 _ 2{x - y)Tr^Y 



^ (6x^ - 3x2 _ _ 37^^2 _ ,^(6x3 - 15x2 - 38x - 24)Xy + 



24y3 ^ ^ 12?/2 

7r2 1 

+ -^(2x^ - 9x - 8) (2x3 - 9x2 _ ^^2^; _ 2)Y^ 

8?/^ 8xy 

+ «7r — ( (3 + 6x - x2)xLi3(-y) + (x^ - 2x2 - x - l)Li3(-x) + 



2y3 



+ 2(x - y)xyLi2(-x) + ^(3x2 + 3x + 2)XLi2(-x) + 

+ ^(x - i)(3x + 2)x^ + ^{x^ + lox + 5);ty2 _ 

7r2 \ 
- y(3x2 + lOx + 5)xy + 3C3ic2 j - 



^ (^xY^ + (x - 1)X^Y + TT^xX^ + 



% V 

+ ^ ( -(6x2 _ 4^ _ i)Li2(-x) - -(39x2 _^ 26x + 4)^2 + 
Ay^ \ 6 

1 7r2 
+ -(11x2 + 20x + 1)XY + -^(2x2 - 8x - 1) 

- -^y2 _ —L(6x3 + 4x2 + 25x + 24)X - 
24y 12y3 

1,0 o , 26X + 23' 

r(6x'^ + 25x2 + 18x + 6)y + 



12xy2 6y2 



(A.67) 



-55- 



J; 



[3] _ X{X - 1) 



Lu{-x) + Lu{-y) - Li4 l^-^j - (Liai-x) - Cs)^ + 
+ (Lisi-x) - C3)Y + yLi2(-x) - ^X^ + lxy3 _ 



24 



6 



12 



120 



+ 



+ 



+ 



3y2 



x(34 + llx - 20x'^)X^ + 3(3x + 2){x^ -x- 3)X'^Y ) + 



+ 9x + 2)(Li3(-a;) - XU2{-x)) 



- (x^ - 6x - l)(Li3(-y) + YU2{-x)) - 

1 

- -(x - 2)(9x + l)Xy2 + — (41x^ - 222; - 15)Y 
8 24 



+ 



TT 



1 



1, 



+ -— (20x + 33)X + -^Cs - + T^(6x^ + 44a;2 + 66x + 25)^2 + 
Soy 6y 9 12y^ 



+ 
+ 



1 



722/2 
1 

2ixy 



{72x^ + 303x^ + 259a; + 22)Xy - (36a;^ + 126x^ + 139z + 73)7r^ 



+ 



(12x^ + 13x2 - 5x + 6)y2 



^ / ^ 19a; - 20,, 
(19x - 17)X — Y + 



72y 



72y 



+ ITT 



2^^3 y 2(8x-^ + llx^ - 2x - 4)XY - (21x'^ + 30x^ + 21x + 14)X^ ) + 
+ ^ (^^^^"^ + + " - |^(x - 25) ) + 



+ ^(153x^ + 230x2 + 19x - 22)X + 

72y3 



1 , q n , 4X — 5 

+ ;:::^^(153x3 + Sllx^ + 16x - 36)Y + 



72xy2 



18y^ 



(A.68) 



x(x — 1) 
6y3 



12(^Li4(-x) + Li4(-y)-Li4(^-^^^ + 

+ 12(y - X){Lh{-x) - Ca) + 27r2Li2(-x) - ^X^ + 



+ 2XY^ - -Y^ - Tr\2X^ + Y^)- — 
2 ^ ^10 



+ 



+ 



3y3 



-(x^ + 15x2 + 15x + 2)(Li3(-x) - XLi2(-x)) - 

- ^(2a;2 - 30x - 27)X^ + ^(Sx^ + 20x2 + 22x + 8)X'^Y - 



- ^a;(x2 - 90x - 87)X + Ca 



+ 



+ 



3y2 



TT 



(x^ + 14x + 1)Y{ Li2(-x) - ^ + 



+ -(7x2 ^ ^ g)^y2 ^ (^2 ^ ;l4^ ^ l)Li3(-y) 



+ 



-56- 



L 



[3] 



X 

X ( \ 



+ 3^(6x3 _^ 23x2 + 32x + Vi)X^ + 



+ -(I2x^ + 35x + 21)XY (27x2 + 81x + 28) + 



36 



+ — (3x2 + 6x-l)y 
4y 



2 x(29x + 9)^ 29x + 



2%2 



+ ivr 



1 

3y' 



3 J Li2(-x) - |(x2 + lOx + 11)X2 + i(2x^ + lOx^ + llx + 4)Xy 



24y 
1 

4' 



-y + 



vr 
24 



(9x^ + 30x2 _^ _^ ^g^J 



3x + 4 2 
24y 



1 



^(x2 - 8x - 5)X ^(x2 - X - 4)y + 

8y^ Sy^ 



7x + 2 
6y2 



in 



(4 + 3x)Xy - —xX^ + -xX - — - —Y' + ^Y + 



Qy\l2'~ ' 12--- ■ 9"-- 2 J 36 

+ (-(7x - 4)X + (7x + 8)1^ - y 

For /i = 5 in eq. ( |2.26| ) and color factor Tr'^1 in eq. ( 2.3lj ): 



54 



(A.69) 



(A.70) 



[1] 



1 

4y 



1-y 



-(20x^ + 17x - 2 + 36x2)Li4(-x) + 
+ 2(18x2 + lOx^ - 1 + 8x)XLi3(-x) + 

- ^(lOx^ + 18x2 + llx + 6)X^ + 
7r2 

+ —(5x3 ^ + + 7)^2 + (5x + 7)C3X + 
3 



X2Li2(-x) - 



vr 



+ ^(1280x3 + 2304x2 + 1105x + 57) 



+ 



+ x(4 + 5x) 



Li4(-y) - Li4 



XLhi-y) 



24 



Y' + 



1. 



1 



+ -Y-'X - lX2y2 _ _y 



vr 



+ 



1 

+ 12 



+ 



1 



6 4 12 

(10x2 + 8x + 2)vr2Li2(-x) - (10x2 + 8x + l)yX(X2 

1 



2vr^ 



(30x^ + 28x - 11) -Li3(-x) + XLi2(-x) + -X^Y 



^(21x + 30x2 _ 25)^3 ^ ^(240x2 + 107x - 205)C3 + 



vr 



(240x2 + 159x - 265)X 



+ 



+ 



1 vr2 

(721x2 + 1449x + 368)X2 - — (lOx - 3) 



288y2 



+ 



-57- 



+ ITT 



5 2y 



^ ( (lOx^ + 18x^ + 8x- l)Li3(-x) + ^—^XU2i-x) 

- ^{lOx^ + 18x^ + llx + 6)X^ - 

- —{5x^ -1 + Ax + 9x'^)X - (x - y)C3 
- x(4 + 5x)Li3(-y) - ^(lOx^ + 8x + 1)X'^Y + 

+ ^ (^(30x2 + 28x - ll)(Li2(-x) + XY) - 

1 vr^ 

- -(SOx^ + 21x - 25)^2 - ^(120x2 + 83x - 73) 1 - ^Cs + 

-(7x3 _ g^2 _ -^^^ _ 5)Li4(-x) - x(7x2 - 6x - 14)Li4(-y) - 



- {bx'' + 12x^ - 4)Li4 (^-- j - 2(3x^ - 3x - 2)yyLi3(-y) + 

+ (Gx^ + 12x2 + X - 4)yLi3(-x) + 

+ (x^ - 12x2 - 7x - l)XLi3(-x) + {x^ + 5x + 4)XLi3(-y) + 
+ ^(6x2 + 3x + l)X'^U2{-x) - (3x2 + + 2)y2Li2(-x) + 

+ y(5x3 + 12x2 - 18x - 14)Li2(-x) - 4yLi2(-x)Xy - 

- ^(llx^ + 8x + 24x2 - 2)Y^ + ^x(2 - x2)X^ - 

- ^^(a;2 - 2x - 2)x2y2 + ^(x^ + 6x2 + 53. + i)^^^ + 

+ |-(llx - 10)Y^X + (4 + 5x)FC3 - ^(9 + lla;)^C3 + 

7r2 7r2 

+ — (x^ - 9x2 _ -^4^ ^ -^2)^2 + _(a.3 _ 3^2 ^ -^2x + 14)y2 _ 

7r2 

- -g-(2a;3 - 9x2 ^ ^^^^ ^ I2)yx + 

+ (488x3 - 1152x2 - 987x - 343) 

1440 ^ ^ 



(A.71) 



1 

Ay 



+ 

2 



x(7x^ + 16x + 21)(Li3(-y) + yLi2(-x)) + 

+ x(7x2 + 2x + l)(Li3(-x) - XLi2(-x)) - 

- —(3x2 ^ ;^g^ ^ 29)X3 - ^(6x2 + 15x + 1Q)Y^ + 

18 3 

2x 1 
+ — (4 + 9x)yx2 + -(3x3 _ ;l2^2 ^ 2lx - 40)Xy2 _ 

O 

1 7r2 

- ^(360x2 + 133x - 443)C3 - — (x^ + x2 + 26x - IQ)Y 



-58- 



a 



[1] 



- ^(96x^ + 72x'^ - 541x - 477)X 



+ 



1 

144^2' 



+ ^(99x2 + 73 + 72x)XY + j^i^^x^ + 81a;^ + 173x + 27)Y'^ 



TT 



288y 



(198x3 + 108x2 - 145x - 455) + 



1 



1 19139 



+ (81x^ + 647x + 782)X - -(1 - x)F + -x + 

216y^ ^ 8^ ' 8 2592 

+ ( (''^^ 6x)Li3(-x) - 2x(3x + l)yLi2(-x) - 

- x(5 - 7x^)Li3{-y) + (6x2 + 7x + 5)XLi2(-x) 
^(6x2 + 6x + l)y3 - ^ (^2 - 6x - 8)^3 + 



+ 



6 



6 



+ (3x3 _ ^ ^ l)xy2 + _(^3 ^ 2x + 1)YX^ - 
6 

- ^(x3 - 12x2 - 5x - 6)X + ^yCs) + 

+ ^ (^Li2(-x)x(7x + 10) + |(3a; - l)(3x + 5)^2 + 

+ -(9x3 ^ 33^2 ^ ^ ;^Q^y2 _ 

6 

- J (9x3 ^ ^2x2 + 13x + 20)Xy + 
7r2 



+ 4^(72x3 + 112x2 + lllx + 79) J - 



1 



2 (27x2 + 12x3 - 16 - 9x)X - 



1 



72yx 



(181x2 + 108x3 + 200x + 27)y + 



2y[24^ 



(x2 - 2) -2ALu{-y) + 24Li4 



1 647 
y " 216. 

47r2Li2(-x) + 



+ -(7 - 2x2)Li4(-x) + (Li3(-x) - C3)y - 

- y{x - l)XU3{-y) + (x2 - 3)XLi3(-x) + jLi2(-x)x2 - 



+ 



(A.72) 



+ — (1 - x)x2y(3y + 2X) + — x(x - i)yx 

24 6 



-59- 



^^{x + 3)XC3 + 7^(32^' - 48x + 11) + 
+ ^(aj + 3) f-SLial-y) - WU2{-x) - ^XY^ + ttV ) + 



+ 



1 

4y 



32y- 



x(7 + ^2 + 2x)(Li3(-x) - XLi2(-x)) - 

- + 6x - 1)^3 + x(x - 2)YX^ + 

1 vr^ 
+ -(4x2 + 7a; + 15)C3 - ^(2x=^ + 6x'^ + 7x- 3)X 

:{2x^ + Wx'^ + 21 + 21x)X^ + ^{x + S)XY - 



1 vr^ 

(x^ + Sx^ + 3x - 3)y2 (6x2 _^ _ + 

16x 96 

3 , , 1 , o ^ x 255 



+ in 



^ ' x(2 - x2)Li3(-x) + |Li2(-x)X - (2x + 3)C3 



2y 



+ 



- + !^x(2 - x^)X^ + J(l - ^) (2Li3(-y) + 1^X2) + 

x(x - 2) 



2y 



U,i-x) + XY-\x^-\)-\.^- 



X ,^ 3x-l),, 1 25 

— -^(7x + 15)X ^ '-Y - — 

16y2^ ^ 8x 2y 32 



(A.73) 



D^^^ = x{x-y) 



-Lu{-x) - hu{-y) + Li4 (^"^ j + 

+ (Li3(-x) + Us{-y))X - yLi2(-x) - ^Y^X + 



1 



12 



24 



TT 

12' 



1 

24 



X^ + 



1, 



+ ^yx^ - —Yx + — TT' 



6 



45 



|Li3(-x) - ^xLi3(-x) - ^(1 - 6x)YX^ - 

6x)XLi2(-x) + 3^(24x2 + 15x - U)X^ 
6 72y 



TT 



65 



17 



TT 



(48x^ + 21x - 47)X + — Csx + — (3^ + 7t;;(4 + 9x) - 



24 



24 



27 



144y 

^(70x2 + Mix + 35)X2 + ^(163x - 53)X + ^ + 



864 



+ iTT 



x(x - y) ({Us{-x) + Li3(-y)) + ^FX^ -^x'^-^X]- 
- ^(1 - 6x)Li2(-x) + 2^(24x2 + 15x - U)X^ + 



-60- 



+ i(6x-l)Xy-^(16x-3)- 

--l^(141x + 70x2 + 35)X-i-l^ 
36y^ y 216 



(A.74) 



x{x — y) 



2yLii{-x) + Ui{-y) - Li4 [-^ 
- 2X(Li3(-x) + Li3(-y)) + yLi2(-x) + U^^'X 

12 3 2 12 6 



TT 



8 



+ -vr^yX - — vr^ 



+ 



3y 



3 

„3 (\J2 



+ 



6 45 

(x'' - - 9x + 2)(Li3(-x) - XU2{-x)) + 



+ ^(8a;^ - 24x - ig)^^ + ^(4x^ + ISx^ + Ibx + Q)XY'^ - 

1 vr^ 
- -(8x^ - 21x2 _ + 8)yx2 - ^(Sx^ + 33x2 _^ 2lx + 14)y + 

vr^ 1 
+ ^(8a;^ + 120x2 + 73x - 27)X + —(288x2 + 287x - 25)C3 

^(Li3(-y) + yLi2(-x)) + ^(3x + 2)y3 _ 

^ (18x^ + 24x2 - 79x - 139)^2 - —(9x2 _ -^g^ ^ -^Ny2 _ 

36 



72y2 
1 

144y 
1 



8x(9x^ - 3x - 7)XY - (36x^ + 39x - 33 + 68x2)7r2 



216,(191 + 335x - 180x2)X + i(5x " 4)^ - || + 



+ in 



|(x - y) - 3yx2 - 6Li3(-x) - 6Li3(-y) + tt^x"^ + 

+ — f (12x2 ^ -^2x - l)Li2(-x) + -(4x2 - 39x - 32)^2 - 
3y V 8 

- i(4x2 + x + l){y- l)y2 - i^(4x2 - 39x - 41)Xy + 
8 4 



+ —(24x3 - 42x2 - 45x + 17) ) 

48 ^ 7 36y 

1 3 371 

— (6x2 j^iQ + — + 

18y^ ' 2y 216 



2(12x2 - 59x - 125)X 



G 



F, 



[1] 



4 1 



G 



[1] 



(A.75) 

(A.76) 
(A.77) 



For /i = 5 in eq. ( |2.26| ) and color factor Tr'2] in eq. ( 2.31| ): 



A 



[2] _ 1 



1 

48' 



(1 - y)(l - xy) 6X2Li2(-x) - 12XyLi2(-x) + 



- 61 - 



{3x-^ + 9x^ + 9a; + 2) -U4 — Li2(-x 



y 



24 



+ 



- ^(3a;3 + 9x^ + 9x + 5)X^Y - -^X{x^ + Sx^ + 3x + 4)(y^ - tt^X) - 

2 c 4 

- - -C3X - Cal^ - ^(4x3 + 12x2 + 12x - 7) 
+ y(^-n(Li3(-x)+Li3(-y)) + 

+ ;5^(6x2 + Ihx + 44)(Li3(-x) + Li3(-y) + {Y - X)U2{-x)) - 



^(23x^ + 60x2 _^ 42x + 50)X^ - ^(2x^ - ISx^ - 63x - 32)X2y - 
- ^(21a:^ + 78x2 ^ -^q^^ ^ 82)Xy2 + 

+ ]^(86x2 + 282x2 + 327x + 208)y2 + 
7r2 

+ — (92x^ + 288x2 + 288x + 335)X + 
288 

J2 3 
+ 777(86x3 _^ 276x2 + 312x + 267)y + —(4x2 + lOx - 13)C3 + 



16 



+ 



288y _ 



(475x3 + 798x2 ^ ^gg^ ^ 368)^2 - 

- 2(475x3 + 1062x2 + 846x + 280)Xy + 
+ (475x3 + 1326x2 + 1506x + 1160)y2 + 

+ 7r2 (475x3 + 982x2 + 622x - 528) 



1 

+ - 

y 



+ m 



-^(950x2 + 1433x + 501)X 
288 



1 

" 864 
9 

4y 
11 
36y 



+ 



36077 475 



(2850x2 + 5091x - 32)y 
C3 - 11 ^^77^^(1 - xy) [{X - y)2 + ;r2) + ^7r2 + 



3456y 288 
(1-y),,. ,.,0 o\ 11 



12y 



+ ^((6x2 + 15x-2)(y-X) + 22y)+|| + 



288y 
2327 11 
86%^ 12 



(A.78) 



[2] 



1 

2yL 



1 



X 



--(y - l)(4x2 - 5x - 5)Li4|^--J - (12x + 13)xLi4(-y) - 
- (12x2 + 9x + 4)Li4(-x) + (x3 + 9x2 + x - i)xLi3(-x) - 



(x3 - 9x2 

(x3 + 3x2 



14x)yLi3(-y) + (x - 5)y2xLi3(-y) - 
4x)yLi3(-x) - 7(9x2 ^ ^ _ 2)x2Li2(-x) 
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TT^ 1 

- — (4x^ + _ _ 26)Li2(-a;) + ^{y - l)(3x + 5)XYU2{-x) + 

+ ^(3x2 + 3x + 2)y2Li2(-a;) + ^(x^ - 3,t - 3)X^ - 

-l{x^ + 3x2 + 2x - l)X=^y + -^(6x3 - 3x2 _ _ 22)X^Y^ - 
6 8 

- ^(lOx^ - 21x2 - 47x - 28)XY^ + 

1 7r2 
+ 24(5^^ + - 2x - 4)y^ - y (x^ + x2 - 7x - 4)Xy + 

7r2 1 
+ — (2x-^ + 6x2 ^ 3^ _ ^2)x2 + -(I2x + 11)C3^ + 

7r2 

+ 24(43;^ - 3x2 _ _ 34^y2 _ ^g^ ^ 5^^^y ^ 

+ (68x2 + 1404x2 + 744x + 343) ^ 

1440 



1 

+ — 
iy 



-x(3x2 + 2x - 4)(Li3(-x) - XLi2(-x)) - 

- x(3x2 + 2x + 8)(Li3(-y) + yLi2(-x)) - 

- :^(41x2 + 24x + 15)X^ + ^{I9x^ - 6x - 4)X2y + 
18 6 

+ l{3x^ + 36x2 + + 40)Xy2 - -^(43x2 + 78x2 ^ ;l20x + 60)y2 
6 18 

7r2 

- —(128x3 + 72x2 + 552x + 477)X - 



7r2 



1 

+ 7% 



^ ^ (25x3 + 54x2 - 75x + 60)y + (144x2 + 36x - 443)C3 
18 36 

^(359x3 + 186x2 + 186x + 144)^2 - 



+ 



(359x3 ^ 3722;2 + 265x + 144)Xy + 

+ 



+ ^{nSx^ + 784x2 + 382x - 167) 

+ (359x^ + 558x3 ^ 344^2 ^ gQg^ _ 27)y2 ^ 

144xy 

+ (915x2 + 402x - 782)X - 
216y 

1,9 , 14135 139 



(6x^ + 5x + 4)Li3(-x) + -(6x + 5)Li3(-y) + 

y ^ 

+ (3x + l)xyLi2(-x) - (3x2 + 3^ + 2)XLi2(-x) + 

+ lx(x2 - 2)x2y + gy3 + + 5)X3 - 

1 7r2 
- -(x3 - 6x2 _ 3^ + 2)Xy2 -—(x^ + 12x2 + 5x + 6)X + 



-63- 



+ - 12^2 - 16x - 8)y + ^Csj + 

+ ^ ^-72a;Li2(-a;) - x{31x'^ + 36x + 7)X'^ + 

+ 2(31x^ + 36x2 ^ 29x + 40)XY - {31x^ + 36x'^ + 87x + 20)Y^ 

2 \ 

- ^(124x^ + 128x2 + 116x + 79) j - t^C^^^x + 79)X + 

1 571 3 1 ~ 

+ ^;^(186a;^ + 79^2 + 164x - 27)y - — — + — 
72xy^ ^ 108y 24 



(A.79) 



1 

2y 



-|(5x2 + 12x + 2)Li4(^-^) - 

- x(6x2 + 6x + l)(Li4(-x) + Li4(-y)) + 
+ X{2x^U3{-x) + 2x^U3{-y) + SCs) + 
+ (4x2 + 6x + l)xy(Li3(-x) + Lhi-y)) - 



- y)2Li2(-x) - ^(9x2 + 12x + 2)Y^ - 



X' 



—X'^i-SXY - 9^2 + 2X2) _^ 

24 ^ ' 

2 

+ ^(5x2 + 12x + 2)xLi2(-x) + ^(15x2 + 24x + A)XY^ + 



12 

+ Il^'(2x2-i2xy + 3y2) + ^ 



(116x^ + 90x2 _^ -^5^ _ -^^^ 



+ 



1 

+ 7^ 



3x(2x2 - 2x + l)(Li3(-x) - XLi2(-x)) + 
+ x(2x2 + 6x + 5) (3Li2(-x)y + 3Li3(-y) - 7r2y) + 
+ ^(6x2 - 2x + 1)X^ - ^(4x2 - 2x + l)xX2y + 

+ 2(3x2 + 4x + 4)xXy2 + l(22x + 13)xy^ - 



I (3x2 ^ ^ _ 3^^2jj^ _ ^^Q^2 ^ -^9^ ^ -^5^^^ 



+ 



+ 



16y 



-(11x2 + 22x - 6)^2 - x(llx2 + lOx + 4)Xy + 



+ y(33x^- 26x2 + 29) 



+ 



+ —(46x2 - 17x + 3)X - —(46x2 - 5x - 12)y + 
32?/ 32y 



+ 



+ ITT 



32xy 
1 /x 



(llx'' - 2x3 - 10x2 - 8x - 6)y2 + 



187 17 

-7:7: + 



128y 32 



yV2 



(6x2 ^ ^ l)(Li3(-x) + Li3(-y)) - 
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TT 



- tt;^ + -ir^ Y - —x^X + -Cs - 



12 



12 



X 



- ^(5x + l)Yi3XY -Y^- TT^) 



1 /3 



, ^ ,x - y)xLi2(-a;) + — (Qx'^ - 2x^ + a; - 3) + 
y V2 24 

+ J (3^2 - 2x + 1)X{X - 2Y)] - J(3x + 7)^^ + 



+ 



X 1 27 3 

8y^ ^ 8x^ ' 32y 8 



(A.80) 



4^1 = J 



-Li3(-x) - Li3(-y) + XLi2(-x) - yLi2(-x) + 

+ ^(5x^ + 15x2 + 15x + 14)X=^ -liy- i)ixy - 1)X^Y + 
24 8 

+ ^iy- l)(xy - l)y3 - -^(3x3 + 9x2 ^ ^ io)xy2 + 
z4 8 

+ ^(lOx^ + 30x2 ^ ^ 37^^ ^ 

+ ;^(x3 + 3x2 + 3x + 3)7r2F-^C3 + 



24 ' 



+ 



72y 



-(37x^ + 60x2 ^ -^5^ ^ 35^^2 ^ 2(37x^ + 72x2 + 45x + 9)XY 



(37x3 ^ 34^2 ^ ^ 7j)y2 _ !L(iiia;3 ^ 210x2 + 120x - 191) 



1 



-(222x^ + 294x - 53)X + ^(37x^ + 55x + 28)y + 



1307 37 



+ ITT 



216y^ " ' ' 36y 

1-2//-, ^..-.l,^ ^^x2 , _2 



864y 72 



^ + 



6y 



-(l-xy)((X-y)2 + r 



144y 



7r^ + 



1 



18y 



+ — (6x2 ^ _ y^^^-^ _ ^g^2 ^ ^ 3^^y ^ 



185 



1 



) 216y 6 



(A.81) 



= x(x — y) 



2Li4(-x) + 2Li4(-y) - 21A^ 



TT 



1 



- 2X(Li3(-x) + Us{-y)) + yLi2(-x) + 

+ ^x^ - ^x^y - ^x2y2 + ^XY^ - 

(x^ - 9x2 _ 9^ + 2)(Li3(-x) - XLi2(-x)) + 

+ ^(7x3 + 9x2 ^ -^5^ ^ g^y3 ^ A2;(l3x2 - 24x - 2A)X^ - 

- l{9x^ - 27x2 - 25x + 8)x2y - -^(y - l)(5x2 + 8x + 2)Xy2 + 



-65- 



+ — (6x^ + AOx^ + 48x + 9)X + 
16 

7^2 1 

+ — (y - + 13x + 1)Y + —(288x2 + 288x - 23)C3 



|-(Li3(-y) + YU2{-x)) + ^(x^^ - 11x2 + 57x + 123)X2 + 



4(x2 + 12x + 31)xXy - 2(x^ + 36x2 _^ _^ 3g^y2 



144y 



(2x3 - 60x2 _ _ 33^)^2 



— (Ilx2 + 8x-36)y + 
36y 



1 , o , 3401 19 



+ ivr 



x(x - y) ( -2(Li3(-x) + Li3(-y)) + ^X^ -X^Y + ^X] + 



+ 



- ( ^(12x2 + 12x - l)Li2(-x) + 7^(x - 5)(8x + 7)X^ - 
y \3 24 

- ^(SX^ - 33x2 _ 3^^ _^ g^)^y _^ 

1 7r2 

+ —(8x3 _^ ^5^2 _^ 2lx + 6)y2 + (48x3 - 6x2 ^ - 6x) ] - 

24 144 ' 



X 2 - 1 



2-(12a^ -7x-46)X 
18y^ 

551 2 
216y 3 



18y 



(12x2 + 17x + 19)y + 



[2] 



4 ' 



[2] 



G 



[2] 



(A.82) 

(A.83) 
(A.84) 



For /i = 5 in eq. ( |2.26| ) and color factor Tr^^l in eq. (|2.31| ) 



rrm _ 1 



-(33x3 _^ gg^2 _^ _ 5)(Li^(_^) + Li4(-y)) - 

- 21Li4(-y) + -(44x3 _^ g3^2 _ 7^ _ io)Li4( -- ) + 

2 V y 

+ (28x3 + 51x2 + 22x - 2)XLi3(-x) + 
+ (24x3 + 39x2 _^ 5^ _ i)XLi3(-y) + 
+ (3x3 _^ g^2 _^ _^ l)yLi3(-x) + 
+ (7x3 _^ 21x2 + 29x + 13)yLi3(-7/) - 

- ^(x + 3)(2x2 + 3x + 4)x2Li2(-x) - 

- ^(2x3 + 3x2 _^ ^-^^ _ 6)xyLi2(-x) + 
+ ^(6x3 + 15x2 + 27x + 10)y2Li2(-x) - 
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2 

- ^(44a;^ + 6Sx^ - Six - 22)Li2(-x) - ^(26a;2 + 33x + 11)X'^ + 

+ ^(24a;-'^ + 27^2 - 6x - 22)X^Y + i(50x^ + 105x^ + 29x + 26)X'^Y^ - 
6 8 

- ^(30a;3 + 39^2 - 592; - 36)XY^ + ^(Ha;^ + ISx^ + 2x- 1)^^ + 
+ —{Ux^ + 27^2 + 26x + 18)X2 - —{52x^ + 93x^ + 33x + 10)XY - 

- — y(50a;2 + 43x - 38)^^ + (g^; + 1)^3^ - (9a; + 2)C3l^ + 



— (985a;^ + 1875^2 + 900a; + 184) 
360 



1 

^ 12y 



+ 

3 «n^2 



(3\x^ - 60a;^ - 20x + 124)(Li3(-a;) - XLi2(-a;)) + 

+ {3\x^ + 114a;2 + 119x + 62)(Li3(-y) + yLi2(-a;)) + 
+ ^(40x2 - 63x - 35)X-'^ - \{3x^ - 30x - 22)y^ - 

- (36x3 _ 242^2 _ -^4^ _^ \9)X^Y + 

+ (32x3 + 102x2 + lOlx + 41)Xy2 + 
+ Y^(98a;3 + 192x2 - 275x - 780)X - 
7r2 

- Y^(112x3 + 588x2 _^ 257x + 307)y + (132x2 + 123x - 62)C3 



+ 72^(135a;3 + 24x2 + 268x + 487)^2 + 
1 

^72yL 



6(45x3 + 14x2 + X + 10)Xy - (135x3 - 34x2 - 178x + 250)7r2 



+ 



+ 



+ ^(135x2 + 60x - 172)y2 + (90x - 121)X + ^(45x - 38)y + 



+ m 



^ ( (31x3 + 60x2 ^ ^ 12)Li3(-y) + 



2y 



+ (31x3 + 60x2 _^ _ i)Li3(_2;) _ 

- (2x3 + 6x2 ^ ^ 3)XLi2(-x) + 
+ 2(l-y)(2-xy)yLi2(-x)- 

- -(14x3 _^ 24x2 + 13x + 4)X3 + 
3 

+ ^(26x3 + 42x2 + 9x - ll)x2y + 
+ ^(5x3 + 18x2 + 22x + 22)Xy2 - 

- ^(3x3 + 12x2 + 9x + 2)y3 - (5x - I)y27r2x - 

- ^(1 - y)(7x2 + lOx + 13)7r2y - C3] + 



-67- 



+ - ( :j^(174x2 + 139x - 62)Li2(-x) + 

+ -^(39x^ - 138a;^ - 62x + 86)X'^ - 
- J(13x^ - 46x2 - 37ar + 6)XY + 
+ J (13x3 + 12x^ - 7x - 8)y2 + 



+ ^(234x3 - 180x^ - 136x + 9) 
: (153x2 - 223x2 - 454x + 30)X 



36y' 
1 



y(153x^ - 115x - 202) + 

36y 



211 
108y 



(A.85) 



r[3] _ 



3 1 

--(x^ - l)xLi4(-y) - -(3x2 _|_ ^ _ i)L^(_a.) _ 



- ^(16x2 + 9x2 + 3^ + 2)Li4^-^^ - ^(2x2 + 3x2 + 2x + 2)XLi3(-x) - 

- ^(6x2 + 3x2 + X + i^xUsi-y) + ^(3x2 + 3x2 ^ + l)yLi3(-x) + 



3x2-7x-4)x2Li2(-x) + 



+ ^(7x2 + 3x2 - X + i)yLi3(-y) - ^(2x2 

7r2 

+ ^(16x2 + 9x2 _ 53. _ 2)Li2(-x) + 

+ -^(x - l)(4x + l)xX^ - A(^ - 2/)(32:2 _ 3^ _ 2^x^Y - 
48 12 

- ^(lOx^ + 15x2 + 13x + 6)x2y2 + ^(30x2 + 21x2 + 5x + 8)XY^ - 

1 TT^ 

- 24(42:^ + + 2x + l)y^ - 24(2^^ + 3x + 10)xX2 + 



+ y(4x2 + 3x + 3)xXy - ^(10x2 + 9x2 + 7^ + 6)y2 ^ 
+ ^(x - y)C3^ - ^(2x + l)C3y + ^(25a:2 - 27x + 10)x 
- ^(6x2 - 3x + 2)y2Li2(-x) + ^(2x2 + x + 2)XyLi2(-x) + 



+ 



--(9x2 ^2x- l)(Li3(-y) + Li2(-x)y) - 
- ^(3x2 + 2)x(Li3(-x) - XLi2(-x)) - 



- 1(3x2 _ 3^ + 2)X^ - ^(26x + 3)x2xy2 + 



+ __(22x^_9a; + ll)a;X2y 



(3x^ - 6x - 5)xX + 



+ ^(20x2 + 4x - l)xy - -(16x + 1)C3 



+ 



-68- 



+ ^^y3_|_ ^ ^q^3 ^a„2 a'^^ o^\v2 



12 
1 



24y- 



:{3x'^ - 16x^ - 62x - 37)X^ + 



+ T^(3a;3 - 26x^ - 47a; - 24)Xy - 
2 

^(x=^ - lOx^ - 12x - 8) - -|- (6a; - 41)X + 



8y 



I2y 



+ 



— (x^ - 12x^ - 6x - 2)y2 - i(3a; - 23)y + 
8x 



+ in 



2y 



6 

(x - l)(x2 + X + l)Li3(-x) + (x^ - 2)xLi3(-y) - 

- (2x^ - 2x - l)XU2{-x) - 2(1 - x^)xU2{-x)Y - 

- - yfX-'Y - ^xX - ^x{7x' - 4)y) + 

1 / X 

+ — ( -(2x - 7)xLi2(-x) - g (IQx^ - 18x + 20)^2 - 



--(19x2-12x-3)y2 ^ 



(19x^ - 6x + 6)x + 



+ |(19a;2 - 18x + 19)Xy j + 



+ 7Tr9 + 11^^ + 34x + 24)X + 

47 1 

+ ^— (7x^ + 7x2 + 6)y , 

12xy^ ^ 12y 12 J ' 



x(x - y) 



-Li4(-x) -Li4(-y) + Li4(^-^^ +X(Li3(-x) + Li3(-y)) 



TT 



-Lio(-x) 



—x^ + -x^y + -x2y2 

24 6 4 



1 n 1 



TT 



TT 



12 



12 



45 



+ 



+ 



3y 



(6x2 ^ 5^ _ 2)(Li3(-x) - XLi2(-x)) + 

+ (x-l)(Li3(-y) + yLi2(-x))- 

- (6x2 + 6x - 1)C3 + |^(3a;2 + 57x + 34)X^ - 

- |(3x + 2)y3 - \{2x^ + 42x2 + 37x + 6)x2y + 
6 8 

1 

+ -(x^ + 3x2 ^ _ 2)xy2 + _(4a; _ z){Zx - 5)y + 



vr 



+ ^(3a;^ + 9x2 _^ ^ gg^j^ 



+ 



+ T7r^(3^^ + 13a;2 - 9x - 25)^2 + 
12r 



(36a;^ + 126x2 + 39x - 22)XY - (18x^ + 90x^ + 72x - 73)7r' 



1 

+ TT^i^x^ + 16x2 + 23x - 6)y2 _ -^(36x + 17)X - —(9x + 5)Y + 
24a; ^ ^ 72y^ ' 18^ ^ 



+ 



x(x — y) 



1 



x(x - y) ((Li3(-x) + Li3(-y)) - \x' + ^X^y - 

- ^ ( (6x2 + 4^ _ i)Li2(-x) -\{x^ + 39^2 + 17x - 14)^2 + 

+ ^(.T'"^ + 39x2 _^ ^ 4^^y ~ ^^^^ + ^^^^ + ^'^^ + ^^^^^ ~ 

7r2 \ 
- ^(3x^ + 69x2 ^ + 25) j - 

3_(6x2 + 273x2 + 317x - 22)X - 

72y2 

3_ (6x3 + 195^2 ^ 124x-36)y- — - — 

72xy^ ' 72y 24 J ' 

2^-Li4(-x) -Li4(-y) +Li4(^-^jj - yLi2(-x) + 
+ 2(Li3(-x) + Us{-y))X - - Y^) + 

+ ^XY{X' - y2) + Ix2y2 + _ y)2 _ f.xY + Avr^ 

-(x^ - 9x2 _ 9^ + 2)(Li3(-x) - XLi2(-x)) - 

- ^(5^2 - 24x - 27)X^ + ^(6x3 - 24x2 - 25x + 8)X^Y + 
+ l{y- l)(7x2 + lOx + 4)Xy2 _ 



7r2 

- ■^x{x'^ + 60x + 60)X - (12x2 + 12x - 1)C3 



+ 



+ ^ (hki-y) + Li2(-x)y - y y^ + i(a;2 - 2x - i)y2 + 



+ ^(2x3 + 2x2 - 7x - 13)X2 + —(4x2 - 2x - 5)XY + 
8y^ 8y 

+ ^x(9x + 8) - ^(20x - 9)X - ^(5x - 2)y + 



+ ITT 



x(x - y) (^2(Li3(-x) + Li3(-y)) - ^X^ + X^Y - ^X^ + 
+ - (-^(12x2 + 12x - l)Li2(-x) + 

y \ 3 

+ ^(3x2 - 36x2 _ _^ 4^^y _ ^(3x3 + 12x2 + 12x + 4)y2 
- ^(3x2 - 36x - 38)X2 - ^(3x2 - 4)(3x - A)\ + 



29 7 



+ ^(2x^ - 7x - 21)X + — (2x^ + 7x + 4)y - — 



8y 



24y 



(A.88) 
(A.89) 



B. Auxiliary functions for two-loop scheme shifts 



In this appendix we present auxiliary functions appearing in eqs. (|4.3| ) and ( [4.4D for the 
shift in the two- loop amplitudes under scheme changes. These functions correspond to the 
(5i?-dependent parts of the 0{e) terms in the one-loop amplitude remainders. They are 
given by, 



M- 



M- 



M. 



1 
1 



M. 



M. 



(l),[3].,5^ 



M. 



M. 



(l),[3].,5fl 



M, 



(l),[l]^,'5fl 



N 



1 x\ 1 

4 6 J 4iV ' 



x{2x — 1) X 



12y 



4yN 



0, 
1 



(B.l) 

(B.2) 
(B.3) 



Ni^xX"^ + y{l - x)X + y2(2 + 3x) + wx(2X - y(3 + x))j + 
+ 1 (-x^X^ - y(l + 3x)X -y\2 + x)- 



(l),[2]e,5fl ^ ^( 

4yV 



— i7rx(2xX + y(l — x)) 
-iV(3 + i7r) + 4(l+i7r) ) , 



1 1 \ 1 



4 6x / 4iV ' 



(i),[2].,5« ^^(2-x) 

^ 12xy 4yN 



0, 
1 



N(x^X^ + y{l + 3x)X + y^(3 2x) i7r(2x^X xy(l - x)) 



1 

N 



xX"^ - y{x - 1)X + (2x + + i7r(2xX - xy(x + 3)) 



«'I^]^'^- = l(iV(X-3)-l(X-l)) 
-4^(^^ + ^-^)' 



^iV(^-3)-^(X-l), 



(B.4) 

(B.5) 
(B.6) 

(B.7) 

(B.8) 
(B.9) 

(B.IO) 

(B.ll) 
(B.12) 
(B.13) 



{l),[2]e,5R _ N 



x^(^{X - Yf + vr^j + X(2x - 1) + q2x + 1)Y 



x + 3 



+ 



- 71 - 



+ — x(l + x)X^ - 2x{l + x)XY + Tr^x{l + x) + X{2x + 1) + 



+ x{l + x)y2 + -{2x + 3)Y + - 



y y 



(B.14) 



(l),[3]e,5fl _ 



(B.15) 
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